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ABSTRACT 


This  report  presents  methods  for  predicting  the  performance- 
determining  aerodynamic  characteristics  and  the  stability  derivative*  c 
transport-type  configurations  employing  the  vectored-thrust/mechanical- 
flap  high-lift  concept.  These  methods  are  suitable  for  preliminary 
design.  The /  have  been  automated  in  a  FORTRAN  IV  computer  program,  for 
which  a  users'  manual  is  included  in  the  appendix  of  this  document. 
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SECTION  I 


INTRODUCTION 


1.1  Background 

The  U.  S.  Air  force's  need  for  modernization  of  its  Tactical 
Airlift  capability  led  to  establishment  of  the  Tactical  Airlift  Tech¬ 
nology  Advanced  Development  Program  (TAT-ADP).  This  program  was  designed 
to  contribute  to  the  technology  base  for  development  of  an  Advanced 
Medium  STOL  Transport  (AMST) . 

The  AMST  must  be  capable  of  handling  substantial  payloads  and 
using  airfields  considerably  shorter  than  those  required  by  large  tacti¬ 
cal  transports  now  in  the  Air  Force  inventory.  If  this  short  field 
requirement  is  to  be  met  without  unduly  compromising  aircraft  speed, 
economy,  and  ride  quality,  an  advanced-technology  powered- lift  concept 
will  be  required. 

The  STOL  Tactical  Aircraft  Investigation  (STAI)  is  a  major 
part  of  the  TAT-ADP,  and  comprises  studies  of  the  aerodynamics  and 
flight  control  technology  of  powered-lift  systems  under  cons  uk  rat Lcn  for 
use  on  the  AMST.  Under  the  STOL-TAI,  The  Boeing  Company  was  awarded 
Contract  No.  F33615-7.1-C-1757  by  the  USAF  Flight  Dynamics  Laboratory  to 
conduct  investigations  of  the  technology  of  the  vectored- thrust  and 
internally  blown  jet  flap  powercd-lift  concepts.  These  investigations 
included: 

o  Aerodynamic  analysis  and  wind  tunnel  testing 
o  Configuration  studies 

o  Control  system  design,  analysis,  and  simulation 

1.2  Objective 

The  objective  of  the  work  reported  here  was  to  develop  con¬ 
venient  and  rapid  methods  for  predicting  the  performance-determining 
aerodynamic  characteristics  and  the  stability  derivatives  of  configura¬ 
tions  using  the  vectored  thrust/mechanical  flap  powered  lift  concept. 

The  methods  are  intended  for  preliminary  design  purposes  and  ease  of 
application  has  been  emphasized. 

1.3  State  of  the  Art  Prior  to  the  STAI 

Early  in  the  STAI,  the  available  literature  and  test  data  on 
vectored  thrust  was  surveyed.  It  was  found  that  the  data  base  for  vec¬ 
tored  thrust  interference  effects  on  transport-type  configurations  was 
almost  nonexistent.  Consequently,  the  "State  of  the  Art  Design  Compen¬ 
dium"  compiled  from  the  information  then  available  consisted  only  of 
procedures  for  estimating  power-off  characteristics  and  the  recommenda¬ 
tion  to  correct  for  power  simply  by  direct  vector  addition  of  the  propul¬ 
sive  forces.  That  is,  interference  effects  were  assumed  to  be  zero. 
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To  fill  the  gap  in  the  data  base,  an  extensive  program  of 
testing  was  then  carried  out  in  the  Boeing  V/STOL  Wind  Tunnel.  The 
results  of  that  program  are  reported  in  Volume  IV  of  the  present  series 
of  documents,  and  are  the  basis  for  the  methods  presented  here. 

1.4  Technical  Approach 

Power  effects  are  described  in  this  report  as  the  sum  of 
forces  and  moments  computed  by  direct  vector  addition,  plus  interference 
increments.  The  interference  increments  were  usually  found  to  be  best 
described  graphically.  That  is,  no  improvement  in  convenience  or  under¬ 
standing  was  apparent  in  attempting  to  reduce  the  curves  to  analytical 
formulae,  except  for  a  general  dependence  of  the  interference  forces  on 
the  square  root  of  the  thrust  coefficient. 

1.5  Scope 

The  scope  of  this  investigation  covers  vectored  thrust /mechani¬ 
cal  flap  high-lift  systems  installed  on  configurations  suitable  for  a 
STOL  tactical  transport.  These  methods  are  intended  to  be  used  in  con¬ 
junction  with  the  USAF  Stability  and  Control  DATCOM  (Reference  1). 

1.6  Document  Organization 

Section  II  presents  methods  for  predicting  performance  deter¬ 
mining  aerodynamic  characteristics  with  power  off,  and  for  estimating 
interference  effects  due  to  vectored  thrust. 

Section  III  presents  procedures  for  computing  stability  and 
control  derivative  corrections  due  to  vectored  thrust. 

The  appendices  provide  a  users'  manual  and  a  listing  of  a 
FORTRAN  IV  computer  program  which  automates  the  procedures  given  in 
Section  II. 
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SECTION  II 


LONGITUDINAL  CHARACTERISTICS 


Aerodynamic  estimation  techniques  are  presented  which  provide 
increments  of  lift,  drag,  and  pitching  moment  for  leading  and  trailing 
edge  devices.  These  increments  are  to  be  added  to  the  clean  airplane 
values  which  may  be  estimated  from  Datcom  or  other  alternate  source. 

2.1  Unpowered  Aerodynamic  Characteristics,  Free  Air 

2.1.1  Lift 

Lift  estimation  below  maximum  lift  has  been  divided  into  lift 
curve  slope  and  flap  lift  increments.  The  effects  of  flap  extension 
(chord  extension)  which  increases  the  wing  area,  and  flap  deflection, 
which  changes  the  wing  camber,  are  treated  separately. 

2. 1.1.1  Lift  Curve  Slope 

There  are  a  number  of  theoretical  or  semi-theoretical  formulae 
which  give  good  agreement  between  the  estimated  and  experimental  lift 
curve  slopes  of  three-dimensional  wings  (Refs.  1,  2,  3,  4).  One  easy-to- 
use  method  is  that  from  Jor.es  and  Cohen  (Ref.  4).  See  sample  problem  for 
additional  definition  of  Sq  and  p.  Page  5. 

-  __  2.TT/A  Sa  _1_ 

(P/bl(AV^  SR  (2'1~1) 


The  modern  high  lift  system  usually  has  trailing  edge  flaps  with 
rearward  displacement  (chord  extension)  and  may  also  include  a  leading 
edge  device  with  forward  displacement.  The  areas  added  by  these  displace¬ 
ments  of  the  eading  and  trailing  edges  must  be  added  to  the  basic  planform 
when  estimating  flaps  down  cLor  t  he  inboard  edge  of  the  flap  is  at 
the  side  of  the  body,  the  added  area  for  flap  extension  will  be  based  on 
the  assumption  that  the  flap  extends  to  the  body  centerline. 

A  comparison  of  estimated  and  test  are  shown  in  Fig.  1 • 
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Data  from  BVWT  097  (Ref  5) 


Figure  1:  Lift  Curve  Slope,  Test  -  Estimate  Comparison 
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SAMPLE  PROBLEM  -  LIFT  CURVE  SLOPE 
STAI  wind  tunnel  model  LE  &  TE  devices  deployed,  15°  sweep. 

SQ  =  Area  ABCDEF  =  8.592  SF 

SRef  -  6.164  SF 
b  =  34.274  in. 

AGross  =  b2/S_  =  5.74  SF 

P  =  ABCDEF  =  100.952  in. 


Calculate  CL^  from  Equation  2.1-1. 

C,  ,  f "ka.vB^  1 

* "  5,1 3 

deg. 

F-eotn  've&T  e>\jucnf  cph  y  eu  s 

2. 1.1. 2  Effect  of  Trailing  Edge  Flap  Deflection 

The  effect  of  pure  (i.e.,  no  area  increase)  trailing  edge  flap 
deflection  is  tr  change  the  zero-lift  angle  (aoL)  without  changing  the 
wing  lift  curve  slope.  The  approach  chosen  here  to  estimate  trailing  edge 
zero-lift  angle  shift  is  due  to  Eldridge  (Ref.  6  and  7). 

Consider  an  infinite  yawed  constant-chord  wing  with  trailing 
edge  flap  defection. 
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& 

It  can  be  shown  that,  referenced  to  the  free  stream  velocity, 

Cf  2.TT  CO&J\  (2.1-2) 

CAh=  StTcX^^CO^A  (2.1-3) 


Therefore: 


corsA 

A-O 


For  flaps  on  tapered  wings,  the  significant  sweep  angle  is  that 
of  the  locus  of  sectional  aerodynamic  centers  for  the  wing,  approximately 
the  quarter  chord  (used  for  C £  ),  and  the  locus  of  sectional  flap  centers 
of  pressure,  approximately  the  half  chord  (used  for  C £$).  If  the  flap 
angle,  6f ,  Is  measured  normal  to  the  hingeline,  then  the  effective  angle 
along  a  chordline  normal  to  the  half  chordline  is 


Se  =  tan  bjf  aos(A  9^  -AhlY] 


(2.1-4) 


AOC 


OL 


2.D 


20  A-O[cos  Acyu.  J 


taV%'ftan^  cc^>(Ac^-/V  H  L)3 


(2.1-6) 


For  a  finite  aspect  ratio  wing,  lifting  surface  theory  shows 
that  the  effective  oig  is  increased  above  the  two-dimensional  value. 
Therefore,  for  wings 


cos(A  c/£  -AmJU  Axe 

(2.1-7) 


Empirical  two-dimensional  data  has  been  correlated  for  single 
and  vane-type  double-slotted  flaps,  Figs.  2  and  3.  Lifting  the  surface 
theory  shows  that  flap  effectiveness  is  affected  by  aspect  ratio.  The  two- 
din.,  isional  test  value  of  a6  can  be  corrected  to  three-dimensional  using  the 
theoretical  results  of  Ref.  8  ,  Fig.  4. 

The  part  span  load  factor  used  in  Equation  2.1-7  may  be  found 
in  Figure  5. 

For  multi-element  clamps,  contributions  of  individual  elements 
add  algebraically  (Fig.  6),  so 


(2.1-8) 
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■e  2:  2-D  Flap  Effectiveness,  Single-Slotted  Flap 


Effectivness,  Double-Slotted  Flap 


Figure  4:  3-D  Effect  on  Flap  Effectiveness 


Element  Flap  Nomenclature 


The  flap  lift  increment,  measured  at  the  angle  for  zero  lift 
of  the  flaps-up  wing  is 


*(cLct  X*  ^°^oi_TEfUxJ) 


(2.1-9) 


Some  flap  lift  will  be  carried  over  onto  the  body.  The  amount 
of  carry-over  will  depend  on  the  wing  position  on  the  body,  the  body  span 
to  wing  3pan  ratio,  and  the  flap  lift  increment.  A  limited  amount  of  data 
for  high  wing  configurations  have  been  correlated  as  shown  in  Fig.  7* 

If  it  is  desired  to  make  the  correction  for  body  carry-over,  flap  lift 
increments  should  be  calculated  assuming  the  flap  ends  at  the  body  side 
and  that  it  extends  to  the  centerline.  The  difference  between  these  values 
is  then  multiplied  by  the  body  carry-over  factor  (k)  from  Fig.  7 
trailing  edge  flap  lift  increment  with  body  carry-over  is 


(2.1-10) 


The  body  carry-over  lift  increment  also  results  in  shifting  the 
angle  for  zero  lift  by 


deacon 


(2.1-11) 


^  ^  OL  -r-p  = 


(2.1-12) 


A  comparison  between  test  data  obtained  from  STAI  wind  tunnel 
testing  and  calculated  data  are  presented  in  Fig.  8. 

SAMPLE  PROBLEM,  TRAILING  EDGE  FLAP  LIFT  INCREMENT 

STAI  Wind  Tunnel  Model  LE  &  LE  Devices  deployed,  T.E.  deflection 
45° / 60° ,  A  c/4  =  15°. 


SG  =  8.952  SF 

SD  =  6.164  SF 
R 

b  =  84.274  in. 
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.4  £  1.2  1.6 

Wing  Flap  Lift  Increment 


Figure  7:  Body  Lift  Carryover 
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Test-Estimate  Comparison 


^Gross  =  5. 74 

p  *  100.952  in. 

Flap  Type  4,  Triple  Slotted,  see  Fig.  6. 


c/C  * 

1.283 

*1 

.323 

Cf2/c'  - 

.091 

6fl  - 

*5.14 

Sf2  = 

15.13 

r‘0B  “ 

=  Trailing  Edge 

'OB  ,/3 

AC'/4  -  15.401 

AC72  =  11.625° 

AHLX  =  7.295° 

AHL2  =  3.525° 

Ct 

a  -  .0988  (calculated  by  method  in  Section  2. 1.1.1) 

Calculate  flap  angles  normal  to  half  chord  line.  Equation  2.1-4. 


6ei  =  tan'1  [tan  43.14  cos  (11.625-7.259)]  =  45.06° 

6e2  *  tan  1  [tan  15.13  cos  (11.625-3.525)]  =  14.98° 

For  forward  flap  section  using  Figure  3,  C'/C'  and  5  read 

ei 

(°‘6)1  =  -•  485 

For  aft  flap  section  using  Fig.  2,  Cr  /C'  and  6 

f  2  e2 

(°i-6 )  2  =  _*25 


From  Fig. 


A  3D 
“*20 


4,  C'f/C 
-  1.03 

1 


and  A  determine 
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From  Fig.  4,  Cf  /C'  and  A  determine 
r2 

-  l-057 
°2D 

From  Fig.  5,  n1B  and  n0R 

Ate  -  .849  -  .183  =  .666 
Since  the  flap  is  the  sum  of  its  parts, 

<A“0L)tE  =  4aOL1  +  4“ol2 
and  aQL  from  Equation  2.1-7 

(Aool)^  C-.485)a.03)(|2|i^il)  (45.06H.666)  -  -14.92 

(AaOL)2  =  (--25)  (1.056)  (14.98H.666)  =  -2.62 

AaOLTE  =*-17.54 

Then  from  Equation  2.1-9 

^lTEflap  =  -(-17. 54)  (.0988)  =  1.73 

body  carry  over  factor  from  Fig.  7,  AC  and  r  (flaps  end 
at  side  of  body)  TE 

K  =  .58 

with  equation  2,1-10 

ACl_f  =  1.73  +  (1.73)  (4fz)  (.53) 


-  1.73  +  .28 
AC.  =  2,01 

lte 


from  test 


CL  =  1.91 
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2.1.1. 3  Effect  of  Leading  Edge  Flap  Deflection 

There  has  been  little  work  done  to  correlate  test  data  on  the 
effect  of  leading  edge  flap  deflection  on  lift  below  ^Lmax*  Since  this 
effect  is  relatively  small  compared  to  trailing  edge  flap  deflection, 
leading  edge  flap  effectiveness  is  taken  to  be  the  potential  flow  value 
given  in  Figure  9. 


2.1. 1.4 


Total  Free  Air  Lift 


The  increments  obtained,  Aciql  and  AC^  and  the  slope  of  the  flaps 
down  lift  curve  may  be  combined  with  flaps  up  estimates  from  Datcom  or  othe 
sources. 


(1) 

ct  flaps  up  from  Datcom  or  other  source 

U  L 

(2) 

A(XOLle 

(3) 

Au0LTE 

(A) 

CXql  flaps  down  vl)  +  (2)  +  (3) 

(5) 

acLle 

(6) 

AC 

lTE 

(7) 

ACLf  =  (5)  +  (6) 

(8) 

CT  flaps  down 

La 

(9) 

Cj  flaps  up  Datcom 

2.1.2 


Maximum  Lift 


Many  attempts  have  been  made  to  develop  methods  for  estimating 
the  maximum  lift  of  an  airplane  with  a  high  lift  system.  No  method  has 
given  consistently  reliable  results.  The  method  given  here  should 
apply  to  the  type  of  configuration  likely  to  be  considered  for  a  STOL 
transport.  Unfortunately,  C^aif  may  vary  widely  from  the  values  calculated 
by  this  method  for  particularconf igurations  with  unusual  arrangements. 

The  approach  taken  divides  the  problem  into  the  Cr  of  the 
clean  wing  plus  increments  due  to  leading  edge  and  trailing  edge  devices. 

For  the  clean  wing,  the  methods  of  Datcom  may  be  used  to  estimate 
.  Next  the  increment  in  maximum  lift  due  to  leading  edge  devices  will 
be  a§aed  to  the  clean  wing,  then  the  trailing  edge  increment  added.  This 
technique  has  been  chosen  because  of  the  availability  of  data  in  this  form. 

It  would  be  more  satisfying  to  add  a  leading  edge  increment  to  the  flaps- 
down  maximum  lift,  since  the  shape  and  optimum  deflection  of  the  leading 
edge  device  is  a  function  of  the  trailing  edge  lift  increment.  However, 
insufficient  data  is  available  to  use  this  approach. 

2.1. 2.1  Leading  Edge  Devices 

The  maximum  lift  increment  due  to  leading  edge  devices  is  a 
function  of  wing  sweep,  device  chord,  shape,  deflection,  2nd  span.  It  is 
assumed  that  care  will  be  taken  in  tailoring  and  fairing  areas  such  as 
intersections  of  nacelle  struts  and  wings,  etc.,  where  relatively  large 
penalties  may  result  from  local  flow  separation  and  interference  effects. 

2 

Correlations  of  ACt  /cos  Ac/4  versus  leading  edge  device 
chord  ratio  are  shown  in  Fig.  Iu  for  conventional  leading  edge  slats 
and  for  shaped  leading  edge  devices  representative  of  current  state-of-the- 
art  variable- camber  Krueger  flaps. 

The  maximum  lift  increment  due  to  the  leading  edge  device  is  then: 

^ C L c°^ A<y4  (2.1-15) 

It  should  be  noted  that  for  this  estimate  the  chord  lengths  are 
measured  normal  to  the  basic  wing  leading  edge  and  that  the  gross  area  is 
the  area  of  the  basic  wing  extended  to  the  body  centerline. 


21 


SAMPLE  PROBLEM,  LEADING  EDGE  MAXIMUM  LIFT 


STAI  Wind  Tunnel  Model,  L.E.  deployed,  15°  sweep. 


.166 

15° 


from  Fig.  10  and  C^^/C  read 


AC, 

^nax 

C0S2Ac/4 


1.14  (shaped  leading  edge) 


with  Equation  2.1-15 


AC,  -  (1.14)COS215.0) 

^naxLE 

=  1.06 


from  test  data 


AC 


WaxLE 


.57 


The  calculated  value  is  too  high  because  the  leading  edge  device 
tested  was  a  compromise  designed  for  several  nacelle  strut  locations  and 
leading  edge  sweep  angles.  A  larger  AC^  for  a  given  configuration  could 
normally  be  achieved  by  tailoring  the  leaiing  edge. 

2. 1.2. 2  Trailing  Edge  Devices 

The  increment  in  Cj^  ^  due  to  deploying  trailing  edge  flaps  is 
caused  by  two  effects;  increaselXarea  due  to  chord  extension,  if  any,  and 
increased  camber.  Assuming  that  the  airfoil  stalls  when  leading  edge 
pressure  distributions  are  similar  for  the  flaps-up  and  -down  cases,  the 
theoretical  maximum  lift  increment  is  related  trailing  edge  flap  lift 
increment  by: 


cairn  be  r 


AC* 


Lac. 


“*KJ7r) 


re 


r*=o 


(2.1-16) 


where 


taken  from  Ref.  9,  is  given  in  Fig.  11. 
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The  maximum  lift  of  a  wing  section  based  on  the  extended  chord 
would  be  almost  unchanged  if  the  trailing  edge  flaps  were  translated  aft 
without  deflection.  When  a  leading  edge  flap  is  used,  however,  increase 
in  wing  chord  would  result  in  a  reduction  in  the  ratio  of  the  leading 
edge  flap  chord  to  the  wing  chord.  This  reduction  in  the  leading  edge 
chord  ratio  reduces  the  increment  in  maximum  lift  due  to  the  leading  edge 
flap  since  this  increment  is  based  on  the  wing  chord  without  extension. 

The  increase  in  maximum  lift  from  a  trailing  edge  chord 
extension  is: 


eVwsnd  flaps  Op 
e.KVe*vston 


lecxdma 

edc^e. 


>+a.3le> 


(2.1-17) 


The  reduction  in  maximum  lift  from  the  reduction  in  leading 
edge  chord  to  wing  chord  ratio  is: 


In  the  foregoing  equations  the  gross  area  is  that  area  of  the 
basic  wing  between  the  outboard  edge  of  the  trailing  flaps  and  the  body 
centerline.  The  trailing  edge  area  is  the  increase  in  the  wing  planform 
area  due  to  chora  extension  with  the  flap  rotated  into  the  plane  of  the 
wing.  The  leading  edge  area  Is  the  increase  in  wing  planform  area  from 
leading  edge  chord  extension  counting  only  the  area  between  the  outboard 
edge  of  the  trailing  edge  flap  and  the  side  of  the  body.  See  sample  problem 
Page  27  for  sketch  defining  areas. 

The  total  Increase  in  maximum  lift  from  the  trailing  edge  flap 
is 


-Vra\V\r\q 


ma 

camba  v~ 


its.  Cl 

'~rr kQLX 

c>\ord 

e^Tensvon 


t-E  clrsord 
r<x  tub 


(2.1-19) 


Figure  12  shows  the  estimated  maximum  lift  coefficient  increment 
correlates  with  test  data  within  -.1. 


SAMPLE  PROBLEM,  TRAILING  EDGE  MAXIMUM  LIFT 
STAI  wind  tunnel  model  L.E.  and  T.E.  devices  deployed,  15°  sweep. 

*  .666  (from  Section  2.1. 1.2) 

TE 

C'f/C'  =  .323 
ag  *  5-74 

AClte  *  .201  from  Section  2. 1.1. 2 
ASte  *  1.104  SF 
^  Gross  «  5.119  SF 
AS'le  -  .624  SF 
Cle/C  a  . 167 
CLE/C'=  .130 
Ac/4  =  15.4° 

*  6. 164  SF 

CLm  =  .98  flaps  up  (»-est  data) 
maxpu 

ACt  =  1.06  from  Section  2.1. 2.1 
^"LE 

from  Fig.  11  and  C'^/C* 

AC^Max 

^  =  -438 

maximum  lift  increment  from  camber,  Equation  2.1-16 

ACI-Max-<-438>(1X7H,<2-01)’1-15 

maximum  lift  increment  from  chord  extension,  Equation  2.1-17. 

AC4ta  '  ('98  +  1'06)<53l9*ifi?624)('666)  '  -27 


From  Fig.  10  at  C^p/C  read  slope  of  curve 

d  ( - ~~) 

COS  A 

c/4  ,  „ 

- - - =  6.9 


27 


Change  in  maximum  lift  increment  for  reduction  in  leading  edge 
chord  ratio,  Equation  2.1-18 

ACLMax  -  (6.9) (COS2  15.) (.130  -  .  167) (---*9 — j ( . 666)  -  -.16 

The  total  increase  in  maximum  lift  from  the  trailing  edge  flap 

ACLMax  “  1-19  +.26  -.16 
=  1.29 

from  test  data 


from  test  data 


^Max 


=  3.12 


The  comparison  between  the  estimate  and  test  data  show  a 
fortunate  combination  in  the  estimated  data.  The  increment  from  the 
leading  edge  was  low  and  the  trailing  edge  increment  high  resulting  in 
a  better  comparison  with  the  total  from  test  data. 


2. 1.2. 3  Leading  Edge  Boundary  Layer  Control 


The  effectiveness  of  leading  edge  blowing  boundary  layer 
control  is  very  configuration  dependent.  For  example,  a  wing  with  large 
regions  of  separated  flow  near  the  leading  edge  would  show  large  improve¬ 
ment  in  maximum  lilt  with  small  amounts  of  blowing  momentum.  The  correla¬ 
tion  to  be  shown  in  this  section  does  not  include  the  effect  of  BLC  as  a 
cure  for  problem  areas;  e.g.,  separated  flow  in  the  wing/racelle  strut 
intersection  region. 


A  correlation  based  on  unpublished  Boeing  data  is  shown  in 
Fig.  13  for  leading  edge  devices  designed  specifically  for  blowing 
applications.  The  upper  curve  is  based  on  configurations  with  uninterrupted 
leading  edges;  i.e.,  no  wing  mounted  nacelles,  and  represents  a  design  goal 
for  a  well-tailored  configuration  with  wing  mounted  nacelles.  The  lower 
curve  represents  the  level  achieved  with  wing  mounted  nacelles  with  no 
additional  system  tailoring.  An  optimized  leading  edge  device  may  achieve 
thelift  levels  indicated  only  to  fall  below  this  level  when  operated  at 
off  design  conditions.  The  curves  should  yield  reasonable,  achievable, 
levels  but  no  generalized  information  is  available  regarding  best  device 
shape  or  deflection  or  in  what  manner  the  blowing  should  be  distributed  on 
the  wing. 
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Blowing  Boundary  Layer  Control  Effectiveness 


Leading  edge  blowing  boundary  layer  control  may  also  result  in 
some  increase  in  trailing  edge  effectiveness.  This  is  i  result  of  the 
thinner  boundary  layer  that  then  exists  ahead  of  the  trailing  edge  flaps. 
Insufficient  data  exists  to  allow  a  rational  correlation  of  this  effect 
to  be  developed. 

SAMPLE  PROBLEM,  MAXIMUM  LIFT  WITH  L.F..  BOUNDARY  LAYER  CONTROL 


STAI  Wind  Tunnel  Model,  Nacelles  On 


CgLE  -  *06 
Ac/4  =  15.40 
from  Fig.  13  and  Cp 


ACLMax 

C0S'4\c/4 


-  1.3 


Maximum  lift  increment  for  leading  edge  blowing 
ACLMax  “  (1-3) (cos2  15.4) 

=  1.21 

from  test  data 

cLm  “  -29 
Hlax 

This  increment  from  test  data  is  much  too  low,  which  may  be  the 
result  of  off-design  operation  of  the  leading  edge  devices,  i.e.  15°  rather 
than  30°  sweep.  Also,  the  model  had  not  been  tailored,  and  there  were 
grounds  to  believe  that  there  was  trailing  edge  separation  adjacent  to  the 
body.  It  would  be  expected  with  proper  refinement  of  the  model  configura¬ 
tion  the  maximum  lift  increment  from  leading  edge  blowing  would  approach  the 
predicted  levels. 

2.1.3  Drag 

The  approach  will  be  to  divide  the  drag  into  the  clean  airplane 
drag,  the  profile  drag  of  the  leading  and  trailing  edge  devices,  the  induced 
drag,  and  the  pressure  drag  of  the  wing.  Clean  airplane  drag  can  be  found 
uy  conventional  methods. 

2. 1.3.1  Trailing  Edge  Flap  Parasite  Drag 

The  parasite  drag  of  trailing  edge  flaps  is  a  function  of  flap 
type,  area,  and  deflection.  An  empirical  correlation  for  slotted  flaps  is 
given  in  Figure  14. 
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Figure  14:  Parasite  Drag  of  Trailing  Edge  Flaps 
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(2.1-20) 


In  this  correlation  the  flapped  area  is  the  area  forward 
(streamwise)  of  the  trailing  edge  flap  with  the  flaps,  leading  and  trail¬ 
ing  edge,  extended  and  rotated  into  the  plane  cf  the  wing. 

2. 1.3. 2  Leading  Edge  Flap  Parasite  Drag 

The  parasite  drag  of  leading  edge  devices  is  a  function  of 
device  area  and  deflection.  Insufficient  data  is  available  to  establish 
an  optimum  leading  edge  deflection  angle.  However,  unpublished  Boeing  data 
indicates  that  at  the  optimum  angle 


aCDp 


(2.1-21) 


where  the  leading  edge  area  is  the  planform  area  of  the  leading  edge  device 
measured  parallel  to  the  device  chord  plane. 

2. 1.3. 3  Change  in  Induced  Drag  from  Trailing  Edge  Flaps 

Deflecting  trailing  edge  flaps  results  in  a  change  in  load 
distribution  from  that  of  the  clean  wing.  Since  the  clean  wing  is  normally 
designed  to  have  a  load  distribution  close  to  elliptic,  the  loading  due 
to  flaps  will  normally  cause  the  load  distribution  to  depart  from  elliptic, 
resulting  in  an  additional  induced  drag.  A.  D.  Young  (Ref.  10)  gives  this 
drag  for  part-span  flaps  proportional  .a  the  square  of  the  flap  lift 


increment 


0.  =k/^ 

C  \Tf  A 


(2.1-22) 


where  K  is  determined  from  Figure  15. 


More  accurate  estimate  of  the  polar  shape  may  be  determined  by 
methods  such  as  that  in  Ref.  1.  However,  these  methods  require  the  span 
loading  to  be  determined. 

2. 1.3.4  Parasite  Drag  Variation  with  Lift 

Both  the  friction  and  pressure  drag  vary  with  lift.  It  is 
impossible  to  estimate  these  variations  precisely,  yet  some  allowance 
should  be  made  for  them.  The  data  from  a  number  of  wind  tunnel  tests  of 
transport  configurations  with  highly  developed  mechanical  high  lift  systems 
have  been  correlated  to  obtain  the  curve  shown  in  Figure  16.  This  curve 
is  intended  to  give  a  reasonable  preliminary  design  estimate  of  the  parasite 
drag  variation  with  lift  with  both  leading  and  trailing  edge  devices 
deployed . 
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2. 1.3. 5  Induced  Drag 


The  drag  due  to  lift  is  estimated  assuming  elliptic  load 
distribution. 


(2.1-23) 


This  is  used  since  the  drag  increments  estimated  in  the  previous  sections 
are  designed  to  account  for  the  departure  from  an  elliptic  load. 


A  comparison  between  drag  estimated  by  the  methods  described  and 
drag  obtained  from  the  STAI  wind  tunnel  test  program  is  shown  in  Figure  17. 


2. 1.3. 6  Leading  Edge  Boundary  Layer  Control 

The  effects  of  leading  edge  blowing  boundary  layer  control  on 
drag  were  obtained  from  the  STAI  wind  tunnel  test  data.  These  data  indicate 
that 

=  “*5Ca‘-G  (2.3-24) 

SAMPLE  PROBLEM,  FREE  AIR  DRAG 
Sf lapped  -  5.577  SF 
SRef  ■  6.164  SF 
6fi  =  44.9° 

<5f 2  -  15.1° 

(<Vl  =  “-50 

«V2  11  -’25 

Cf^/C'  =  .289  (includes  heading  edge  extension) 

SLE  =  .882  SF 

ClTEf  "  1,73 
Aq  *»  5.74 
A  =  8.00 

"IB  -  *145 
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mi  ;,i , 1 1  . a . —•••-»* 


■ - -  ^..j».ill»iii 


C  ■  3.33  from  Section  2. 1.2. I  and  2. 1.2. 2 
LMax 

C  =  .98  test  data  flaps  retracted 
Max 

Cp  =  .0600  test  data,  flaps  retracted 


“  0 


Calculate  equivalent  flap  angle.  Figure  14. 
6  =  45.14  +  (-^|5)  (15.13)  =  52.94 


Read  from  Figure  14  and  C  *  •  and  ^ 


AC, 


=  .0395 


DP 


TE  C* 


AC 


DT 


TE 


fAC 


=  1.25 


TE  C 


“  .25 


Trailing  Edge  parasite  drag,  Equation  2.1-20 
AC. 


T,  -  (.0395)  <f~~)  (1.0)  «  .0447 
PMIN 


Leading  edge  parasite  drag,  Equation  2.1-21 

;'CD„  -  (.154)  (~f>  '  -0220 

PMIS  6'1M 


From  Figure  15  and  A q  r'1B  and  r'0B  read 


K  =  1.05 
a 

Kf  =  .4 

K  =  (1.05)  (.4)  -  .42 

Change  in  induced  drag  from  t  ailing  edge  flaps,  Equation  2.1-22 


Ac 


(1.73/ 


"Di-  <-42>  •  -°5oo 

at  a  C^  =  2.4  parasite  drag  variation  with  lift,  with  CT 


-  3. 


Max 


D_  =  .0005 
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.2290 


Induced  drag 

c  .  <2-4>2  ■ 

Di  00(8) 

Total  Drag 

CQ  =  .0600  +  .0447  .0220  +  .0500  +  .0005  +  -2290 

“  .4062 

From  Test  Data  at  C  =2.4 

Li 

C’D  =  .3880 

2.1.4  Pitching  Moment 

Deflection  of  leading  and  trailing  edge  devices  affects  the 
tail-off  airplane  pitching  moment  characteristics  by: 

(1)  Moving  the  aerodynamic  center  location  if  chord  extension  is 
involved. 

(2)  Changing  the  pitching  moment  at  zero  lift  because  of  a  change  in 
camber. 

An  additional  effect  which  influences  the  tail-on  pitching 
moment  is  the  change  in  the  downwash  field  behind  the  wing.  In  the  follow¬ 
ing  sections  these  effects  are  examined.  The  methods  for  estimating  the 
change  in  aerodynamic  center  location  and  pitching  moment  at  zero  lift 
are  taken  from  Ref.  7.  Methods  for  predicting  the  effects  of  high  lift 
devices  on  the  downwash  field  behind  the  wing  are  from  Ref.  1  and  11. 

2. 1.4.1  Aerodynamic  Center  Shift  Due  to  Leading  Edge  Devices 

Leading  edge  devices  without  chord  extension  do  not  move  the 
aerodynamic  center  as  long  as  the  flow  remains  attached.  When  chord  exten¬ 
sion  is  present,  the  a.c.  shift  may  be  calculated  by  considering  the 
leading  edge  planform  extension.  The  estimate  of  the  aerodynamic  center 
shift  is  made  relative  to  the  a.c.  location  of  the  basic  trapezoidal  wing. 

An  elliptical  additional  span  load  is  assumed  for  the  trape¬ 
zoidal  wing.  The  part  span  load  of  the  wing  panel  where  the  chord  is 
extended  is  X.  Using  the  Schrenk-Thorpe  span  load  approximation,  this  panel 
load  increases  by  half  the  fractional  area  increase  upon  addition  of 
the  chord  extension  covering  a  small  fraction  of  the  wing  span.  As  the 
chord  extension  tends  to  full  span,  the  panel  load  increment  approaches 
the  fractional  chord  extension. 


The  inner  wing  panel  loads  are  assumed  to  be  centered  at  50 
percent  of  the  panel  span  on  the  local  aerodynamic  center  both  for  the 
original  trapezoidal  wing  and  the  modified  wing.  The  part  of  the  wing 
planform  contained  within  the  body  plan  view  is  treated  in  a  similar 
manner,  letting  the  local  load  move  forward  (or  aft  for  trailing  edge 
devices)  as  dictated  by  the  chord  extension,  but  the  load  on  the  body  is 
held  constant. 

Two  different  equations  have  been  derived,  for  the  leading 
edge  devices  extending  to  the  side  of  the  body,  and  for  outboard  devices 
which  do  not  extend  to  the  body. 

In  the  following  analysis  it  is  assumed  that  the  ba^ic  trape¬ 
zoidal  wing  aerodynamic  center  position  (x  )  is  known  (see  Ref.  1  or  2) 
md  the  value  of  the  load  is  unity,  i.e.  p 

L  *  1.0  (2.1-25) 


Using  Figures  5  and  18  and  taking  moments  about  A  -  A  zives 

i'Y*  ' 


M  -  A  j  Xj  ~ 


(2.1-26) 


In  Equation  2.1-25  x.  is  the  moment  arm  to  the  local  aerodynamic 
center  of  the  trapezoid  where  it  intersects  the  body  (use  Figure  19 
for  correction  to  quarter  chord  location)  and  x„  is  the  moment  arm  to  the 
midspan  of  the  wing  panel  with  leading  edge  devices. 


Assuming  that  M  is  a  linear  function  of  L  and  M  f>  0  at  L  1  0 


leads  to 


3M  _  21 

3L  L. 

Since 

it  follows  from  Equation  2.  -25,  2.1-26,  and  2.1-28  that 
l  .  \  v  .  j'v' 


(2.1-27) 


(2.1-28) 


(2.1-29) 


=  (*a4tvxx 


The  load  of  the  wing  with  leading  edge  devices  extended  is 

*5 

L  »  \  *  aL 


(2.1-30) 


(2.1-31) 


where  y^  —  X2  is  the  load  increase  based  upon  the  Schrenk-Thorpe  span  ]oad 

assumption.  The  area  increase  AS  is  shown  in  Figure  18  and  the  load  factor 
is  obtained  from  Figure  20. 
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Figure  18:  Nomenclature  for  ac  Location  with  Inboard  L.E.  Devices 


Aerodynamic  Centers  Near  Mi 


o  0.2  0.4  0.6  0.8  1.0 

IB 


Figure  20:  Load  Effectiveness  of  Part  Span  Chord  Extensions 
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Again  taking  moments  about  A  -  A  gives: 


(2.1-32) 


In  Equation  2.1-32  x!  is  the  moment  arm  to  the  local  aerodynamic 
center  of  the  wing  panel  with  the  extended  chord  at  the  body  side  and  x^  is 
measured  to  the  extended  chord  wing  local  aerodynamic  center  at  r\ivV\.7- 

Z. 

The  aerodynamic  center  for  the  wing  with  the  leading  edge 
devices  extended  follows  from 


an  _  f  v 

5l  ~  ^ 


/ 


4-N^\l-wu^ 


/ 


•V‘ 


~  nlr\B2. 


Substituting  Equation  2.1-30  into  2.1-33  gives: 


(2.1-33) 


*S\  f  1 

1  X  z  ~  *  2J  (^qc)iTCLp 


X|(X\-X.')+-N 


1 


(2.1-34) 


The  analysis  for  the  outboard  leading  edge  devices  is  very  similar  to  the 
one  employed  above.  A  simplification  here  is  that  the  load  on  the  body 
region  does  not  require  separate  identification. 

Assuming  again  a  unity  load 

L  =  1.0  (2.1-35) 


and  taking  moments  about  A-A  (see  Figure  21)  leads  to 


n 


Hence ,  ». 

KcWP-|u=^^-MC-nC 


(2.1-36) 


(2.1-37) 


and 


(2.1-38) 
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Area  AS 
/  / 


Local  A.  C.  With 
Chord  Extension, 
0.25  Chord 


L  -  l  h  2_ 


(2.1-39) 


Tlie  moment  is 


(2.1-40) 


Using  the  Equation  2.1-28  and  2.1-40  gives  the  a.c.  .location 


frcLC.) 


o.cjlE 


=  SeJfr-V 


\  AO 

Si. 


(2.1-41) 


Substituting  Equation  2.1-38  and  2.1-41  gives 

&S  \ 


(Kac\_  - 

V  I+-U.- 


*2  -  *2J 


1+ 

S>2_  Z- 


(2.1-42) 


and  the  a.c.  shift  due  to  the  leading  edge  devices  becomes  for  the  two 


cases: 


(a)  Leading  edge  devices  extending  to"TlTe'm8Td'e  of  the  body: 


(AXac\ 


AiU\ 


- ^  +  Azjl I  S*)  *2-~  *zj 


(2.1-43) 

-tK«A  „ 


'  ^Z. 


(b)  Outboard  leading  edge  devices: 


(k}Cac)LE  ^ (^trog  _(Xfluc) 


l-t-Xf  £S  *2. 

Sz. 


trctp  (2.1-44) 
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SAMPLE  PROBLEM,  LEADING  EDGE  EXTENSION  ac  SHIFT 


nx  =  .145 

n2  =  1.0 

X  =  35.82  in. 

=  34.76  In. 
X2  =  39.63  in. 
X2'  -  38.38  in. 
AS  =  .882  SF 
S2  =  4.945  SF 


CVTra„  -  57-98 

from  chart  Figure 


20 


Hs  =  .985 

from  Figure  5  read 
\1  =  .183 
X2  =  .817 

with  Equation  2.1-43  calculate  shift  on  ac  with  leading  edge  extension 


^ac^LE 


(.183)  (34.76  -  35.82)  +  .817  Q1  *  (,985)  (4!945)j  (38,38)  39,63}  4-  37.98 


1  +  (.985)  tfcffP-- 

=  -.194  +  4.488  -I-  37.98 
1.143 


*)  (.817) 


■  36.98  in. 


For  the  change  in  aerodynamic  center,  Equation  2.1-43 

(A  X  )T„  =  36.98  -37.98  =  1.00  in. 
ac  LE 

2. 1.4. 2  Aerodyne  :  Center  Shift  Due  to  Trailing  Edge  Flaps 

Simple  hinged  flaps  do  not  affect  the  aerodynamic  center  sub¬ 
stantially  so  long  as  the  air  flow  remains  attached.  Flaps  with  chord 
extension  move  the  aerodynamic  center  back.  Their  effects  may  be  deter¬ 
mined  by  the  methods  developed  for  the  leading  edge  devices  in  the 
preceding  section.  By  using  the  appropriate  values  from  Figures  19 

46 


and  21  with  y_  from  Figure  20  and  X  from  Figure  5  in  Equation 
2.1-43  and  2.1-44  the  a.c.  shift  due  to  trailing  edge  flaps  extending 
to  the  side  of  the  body  and  outboard  trailing  edge  flaps  may  be  computed, 
respectively. 

SAMPLE  PROBLEM,  TRAILING  EDGE  EXTENSION  ac  SHIFT 

f|^  =  .145 

n2  =  *75 

Xi  -  35.82  in. 

X^  =  37.38  in. 

X2  =  38.22  in. 

X2'  =  39.00  in. 

AS  -  1.104  SF 


2. 1.4. 3  Pitching  Moment  at  Zero  Lift  Due  to  Trailing  Edge  Flap 

is  calculated  by  estimating  the  spanwise  and  chordwise 
position  of  the  center  of  loading  induced  by  the  flap.  AC^^  is  then 
equal  to  the  estimated  flap  lift  increment  times  the  arm  from  the  estimated 
flap  load  center  to  the  flaps  extended  aerodynamic  center. 


The  flap  load  center  'a  estimated  as  follows: 

(a)  As  a  first  approximation,  assume  the  flap  load  center  is  along  the 
wing  half-chord  line. 

(b)  Along  chordwise  cuts  normal  to  the  half-chord  line  evaluate  the 
flap  chord  ratio  c 'f/c»- 

(c)  Determine  the  locus  of  chordwise  flap  load  center  positions  using 
Figure  22. 

(d)  Iterate  if  the  new  flap  center  of  pressure  line  differs  greatly  from 
the  initial  approximation. 

(e)  Locate  the  flap  load  center. 

Finally, 


SAMPLE  PROBLEM,  PITCHING  MOMENT  AT  ZERO  LIFT  DUE  TC  TRAILING  EDGE  FLAP 
Cf/C  =  (constant  %  chord  flap) 


c  =  11.179  in. 
A  *  6.786 
P  =  98.831  in. 

niB  =  *145 

n0B  =  *75 
ACLTe  =  2.01 

Calculate 


Lq  ^  6.786 _ 

C~'  /9aS31\  (6.7861  +  2 
a  \84*>74/ 


.681 
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Figure  22:  Chordwise  Center  of  Load  Due  To  Flaps 
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Front  Figure  22  obtain  chordwise  center  of  pressure 
X 

-  .44,  this  is  near  enough  to  the  original  assuatictn  of  Xcp/  *  .50 

C  c 


that  iteration  will  not  result  in  a  significant  change.  From  Figure  23 
determine  spanwise  center  of  pressure 


.422 


The  center  of  pressure  is  then  located  at  the  intersection  of  the  .44  chord 
line  and  ncp  *  .422.  In  the  model  longitudinal  reference  system 


Xcp  TE 


41.47 


The  (Cm0L)TE  ls  then  calculated 


from  equation  2.1-45 


<Sl>te  -- 1TI79  <«•«-*.») 

-  -.4637 


2. 1.4. 4  Pitching  Moment  at  Zero  Lift  Due  to  Leading  Edge  Devices 

The  pitching  moment  at  zero  lift  increment  due  to  leading  edge 
flaps  is  much  smaller  than  that  due  to  trailing  edge  flaps  so  that  a 
simpler  approach  can  be  adopted 

<2‘1-46) 

Where  X'c/4  is  the  quarter  chord  of  the  mean  aerodynamic  chord  determined 
with  the  leading  edge  extended  in  the  plane  of  the  wing. 

SAMPLE  PROBLEM,  PITCHING  MOMENT  AT  ZERO  LIFT  DUE  TO  LEADING  EDGE 


C.  =»  -.  15 

C  -  11*179  in. 


X' 


C/4 


36.867  in 


X 


ac 


LE 


36.98  in. 


calculate  (AC  )TT, 
m  ,  LE 
ol 


with  equation  2.1-46 
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2.1.4. 5 


m0l  LE  =  jJ-jjI  (36.867  -  36.980) 

=  +  .0013 

Cliange  in  Downwash  Due  to  Leading  Edge  and  Trailing  Edge  High- 
Lift  Devices 


Generally  accepted  methods  for  predicting  the  downwash  variation 
behind  the  wing  due  to  leading  and  trailing  edge  high  lift  devices  are  at 
present  not  available.  However,  qualitative  design  guidelines  based  on  the 
analysis  of  large  amount  of  experimental  data  are  summarized  in  Ref.  11. 
Quantitative  data  for  estimating  the  increment  of  downwash  due  to  trailing 
edge  flap  deflection  have  also  been  obtained.  All  of  the  following 
discussion  is  based  on  Ref.  10. 


Analysis  of  the  air  flow  characteristics  behind  sweptback  wings  shows 
that  before  separation  occurs  the  downwash  remains  unaffected  by  leading 
edge  flaps.  The  increments  of  down  wash  due  to  deflecting  trailing-edge 
flaps  on  wing-body  combinations  are  summarized  in  Figure  24  .  The  ratio 

of  measured  effective  downwash  increment  to  the  factor  aC^ _ was 

found  to  give  satisfactory  correlation  of  the  flap  span  effect  and  is 
shown  in  Figure  24  as  a  function  of  height  of  the  horizontal  tail.  Only 
the  lift  increment  due  to  trailing  edge  flap  deflection  is  used  in  Ref.  10 
indicates  that  leading  edge  devices  have  negligible  effect  on  downwash. 


The  correlation  of 


indicated  in  Figure  24  was 


*  nw  wi.1.  w •*-  k  /  - - v  xuuxuauuu  xu  x  xkuxc  —  ■  was 

found  satisfactory  as  long  as  Ac  was  smaller  than  10°.  When  Ac  was  larger 
than  10° ,  at  low  tail  positions  (close  to  the  wing  wake) ,  the  correlation 
was  not  as  good. 


SAMPLE  PROBLEM,  CHANGE  IN  DOWNWASH  FROM  TRAILING  EDGE  FLAPS 


Zt  =  17.566  in. 
b  =  84.27  in. 


Dig  =  *145 

n0B  =  -75 

ACt  =  2.01 
JTE 
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Figure  24:  Change  in  Down  wash  at  Horizontal  Tail 
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from  Figure  24  read  @  2h/b  -  .417 
he 


hC 


-  15.2 


L/A(n0g  -n1B) 


Ae  =  (15.2)  (2.01)  ,  , 

6  (8)  (.75-  .145) 


2. 1.4. 6  Total  Free  Air  Pitching  Moment 

The  increments  in  zero  lift  pitching  moment  and  aerodynamic  center 
from  extension  are  combined  with  the  flaps  up  data  and  provide  pitching 
moment  as  a  function  of  lift  coefficient. 


OU-  v  ol 

flaps  Ifcadm^  frailmc^ 

up  ex^c^e.  edKje. 


(2.1-47) 


*ac  ~  ^ac  ♦“  ^  Xa.c  +  ^  Kouc.  (2.1-48) 

flaps  learimq  -Vmmrva 

up  ed«je. 


C-m~  ^mOL+' 
flaps 
cAouom 


P 


LC 


etp 


Xqc  "] 

flaps 

douon 


(2.1-49) 


Figure  25  compares  pitching  moment  estimated  with  the  test  data. 

SAMPLE  PROBLEM,  .’ITCHING  MOMENT 


(C 


mo 


flaps 

up 


-.11 


(ACmo>LE 

=  +.0015 

<ACmo>TE 

-  -.4657 

(Xac^  flap 

=  37.98  in 

up 

^ac^LE 

=  -1.00  in 
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+  .820  in. 


(Xac*TE 


X 


37.98  in. 


eg 

from  equation  2.1-47 


(.25  mac) 


(C  )  -  -.11  +  .0015  -  .4657 

mo  flap 

down 

-  -.5742 

from  equation  2.1-48 

(X  )  ,,  =  37.98  -  1.00  +  .820 

ac  f lap 

down 

*  37.80 


Pitching  moment  from  equation  2.1-49  @  a  lift  coefficient  =  2.4 

c„--™+(£?8  37-8° 


.179  11.179 


;)  (2.4) 


-  -.5742  +  .0386  =  -.5356 


From  Test  Data 

C  -  -.630 

m 

2.2  ('round  Effect 

Proximity  to  the  ground  affects  the  wing  aerodynamic  character¬ 
istics  in  three  ways.  There  is  a  reduction  in  dynamic  pressure  at  the 
wing,  a  reduction  in  induced  angle  of  attack,  and  an  induced  camber. 

The  assumption  is  usually  made  (Ref.  1)  that  the  effects  of 
reduced  q  and  induced  camber  are  small  and,  since  they  are  of  opposite 
sign,  can  be  ignored.  While  this  assumption  was  reasonable  prior  to  the 
advent  of  modern  high  lift  systems,  it  is  certainly  not  valid  with  today's 
very  high  lift  ST0L  systems. 

A  very  simple  analysis  has  been  performed  using  a  single  horse¬ 
shoe  vortex  and  its  image  in  the.  ground  plane.  This  will  gi,re  a  theoretical 
estimate  of  the  induced  change  in  angle  of  attack  and  the  reduced  dynamic 
pressure.  The  camber  effect  is  assumed  to  be  small  compared  to  those 
effects  for  ST0L  configurations  with  high  mounted  wings. 

2.2.1  Lift 

To  approximate  an  elliptically  loaded  wing  by  a  single  rectan¬ 
gular  vortex,  the  vortex  span  should  be  Tfb/4.  In  this  analysis  a  single 
horseshoe  vortex  with  span  irb/4  is  used  with  the  induced  velocities 
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averaged  over  the  span.  Consider  the  longitudinal  velocity  (see  Figure  26) 
induced  at  any  point  along  the  wing  span  by  the  image  bound  vortex: 


V  (x'i  -  £^v/CC>2,0v+  C°^ez) 


(2.2-1) 


It  may  be  shown  with  the  assumption  of  the  feme  lift  coefficient  based 
on  the  local  dynamic  pressure  in  free  air  and  in  t round  effect  that  the 
velocity  ratio  is 


Vfa  =  l 


(2.2-2) 


The  ratio  of  lift  coefficients  must  then  be 


BlG£  _ _ * 

I  +  1  + 

tt^A  IL 


\Z1V2 


(2.2-3) 


This  lift  ratio  is  achieved  a«-  a  reduced  angle  of  attack  due 
to  the  induced  velocities  from  the  image  trailing  vortices.  The  change 
in  angle  of  attack  is 


'fa  vf*  -  ^Ojjq 


(2.2-4) 


and  it  may  be  shown  that 


Jr\.  A  ^(i^)  J 


(2.2-5) 


2.2.2  Drag 


The  ratio  of  drag  in  ground  effect  to  that  in  free  air  is 
r-  (  uJr  vu>-x  \l  $ 


'P&e.  _ 


K-Dq  ■'•C'k-* 


(2.2-6) 


C\_P/!s 


(2.2-7) 
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2.2.3 


Pitching  Moment 


The  simple  horseshoe  vortex  approximation  cannot  be  used  to 
find  how  the  center  of  pressure  of  the  wing  changes  from  free  air  to  ground 
effect.  This  would  require  a  more  sophisticated  lifting  surface  analysis. 
As  a  first  approximation  we  will  assume  that  the  location  of  the  center 
of  pressure  does  not  change  in  ground  effect.  Therefore, 


^mCE  __ 


(2.2-8) 


while  this  approach  does  not  have  any  theoretical  justification,  it  does 
correlate  well  with  the  test  data,  see  Figures  27  and  28. 

2.2.4  Downwash 

Using  a  similar  analysis  to  that  for  Cl,  a  and  Cd  it  may  be 
shown  that  the  change  in  downwash  at  the  horizontal  tail  in  ground  effect 
is 


^€GE 


4-  (aH  -zjz)(  4.(2.h--g^  '/s) 


1  fcH- 


i  + 


ft 


pfy2-  vfl-b*- 

L  c*\\ 


J'/a)} 


(2.2-9) 


(ef\D*) 


A  comparison  of  free  air  test  data  corrected  for  ground  effects, 
and  test  data  in  ground  effect  is  shown  in  Figures  27  and  28. 

SAMPLE  PROBLEM,  LONGITUDINAL  CHARACTERISTICS  IN  GROUND  EFFECT 

CL  "  2*° 
h/b  -  .209 

b  -  84.274  in. 

A  *  8.0 

a  =  2.33° 

CD  -  .3410 

Cm  =  -.5273 
m 

Zt  =  17.566  in. 

£t  =  49.171  in. 
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TE  Flip  Span  =  0,743 
TE  Flap  Angle  *  35° 

LE  Flap  Angle  =  70° 

C  /i^  ~  0.06 

Data  from  BVWT099  (Ref.  5) 


For  lift  if  ?.ouBd  effect,  equation  2.2-3 


^GE 


For  angle  of  attack  in  ground  effect  2.2-5 

•  ^  ^5  -  i.qo 

*  .93 


For  drag  in  ground  effect,  equation  2.2-6 


4  -2J&2L 


Pitching  moment  in  ground  effect,  equation  2.2-8 


Change  in  downwaah  in  ground  effect,  equation  2.2-9 


6£ne  — - 


■*"  ^.nf^z-Xn^O  -naif- 


S.  -,wcn  EAD 

=.  -  fe3H-° 
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2.3 


Vectored  Thrust 


This  section  contains  formulae  for  longitudinal  force  and 
moment  coefficients  incorporating  thrust  effects  and  a  discussion  of 
thrust  interference  effects  on  these  coefficients.  The  longitudinal 
force  and  moment  coefficients  are  presented  below. 


Cl  ~  -Cj  w>(«vO  (2‘3-1> 

ore 

~  C  °Ffc*«e.e  +'^t>)K{T  ~  Ct  COS(^  +  <*'}  +  Co^  (2.3-2) 

OFF 

Cm  -  S»n<r+-  co->d)  (2.3-3) 

ore 

+  Co«am(J^- 


The  interference  effects  presented  were  obtained  from  the  STAI 
wind  tunnel  test  BVWT  099.  These  effects  are  the  differences  between  the 
power-on  and  power-off  test  data  with  the  appropriate  thrust  component 
removed  from  the  power-on  data.  The  interference  corrections  are  shown 
as  functions  of  thrust  vector  angle,  angle  of  attack,  nozzle  longitudinal 
location  and  nozzle  grrss  thrust  coefficient. 

The  vertical  and  spanwise  location  effects  are  apparently 
negligible,  although  the  available  data  was  limited.  Spanwise  locations 
tested  were  from  27%  to  60%  of  wing  serai-span.  The  nacelle  centerline 
heights  tested  were  h/c  ■  .371  and  .453  below  chord  plane.  These  vari¬ 
ables  are  not  included  in  the  estimating  procedure. 

The  vectored  thrust  interference  data  were  analyzed  to  general¬ 
ize  the  data  with  sufficient  accuracy  for  preliminary  design  purposes. 

The  methods  will  provide  good  results  for  configurations  having  reasonably 
high  aspect  ratios  and  engines  located  under  the  winji*  since  the  data  base 
for  their  derivation  was  so  restricted.  Application  to  other  arrangements 
is  subject  to  considerable  uncertainty.  Figure  29  shows  the  satisfactory 
agreement  between  measured  forces  and  those  predicted  by  the  present  methods 
which  can  be  expected  when  this  restiiction  is  observed. 

2.3.1  Lift  Interference 

Since  the  chordwise  position  of  the  exit  centerline  of  the  nozzle 
varies  with  vector  angle,  the  d;  ta  had  to  be  crossplotted  to  obtain  all  of  the 
vector  angles  at  the  same  chordwise  position.  The  limited  data  on  spanwise 
location  effects  indicated  that  these  were  minimal. 
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Free  air  lift  interference  due  to  vectored  thrust  may  be  found 
for  x/c  *  .35  and  Cj  «  2.0  from  Figure  30  .  An  increment  for  other 
nozzle  locations  may  be  found  from  Figure  31  .  A  parameter  which  has 
proved  of  some  use  in  correlating  vectored  thrust  and  V/STOL  aerodynamic 
interference  effects  is  the  equivalent  jet  velocity  ratio, 

1/2 

V  is  directly  proportional  to  (1/CT)  .  It  was  found 

C  J 

•  I/? 

interference  correlated  directly  with  Cj  with  sufficient 
accuracy  for  preliminary  design  purposes,  though  it  begins  to  break  down 
at  high  thrust  coefficients  or  angles  of  attack. 

The  lift  interference  for  any  Cj  and  chordwise  nacelle  loca¬ 
tion  is  then 


^r^iNT  L  (2.3-4) 

(f\Ca  (FVG  3f) 


V . 

6  \qjet J 

that  the  lif 


For  this  analysis  the  nacelle  longitudinal  location  is  meas¬ 
ured  from  the  leading  edge  of  the  local  wing  chord  at  the  engine  center- 
line  location,  to  the  center  of  nozzle  exit  plane. 

Symmetric  thrust  conditions  have  been  assumed  for  the  lift 
interference  design  charts  developed.  For  nonsymmetric  thrust  conditions, 
the  charts  developed  may  be  used  assuming  that,  each  wing  operates  independ¬ 
ently  of  the  other.  If  one  wing  has  Cj  *  X  an’  the  other  Cj  *  Y,  the  con¬ 
figuration  will  then  have  a  lift  interference  given  by 


^  =  Z^3+‘2.t^'L\KT®^"'S“  (2.3-5) 


2.3.2  Drag  Interference 

At  a  given  nacelle  location  and  nozzle  vector  angle,  the  free 
air  drag  interference  could  be  correlated  directly  with  the  free  air  lift 
interference.  This  permitted  a  relatively  simple  procedure  to  be  used. 
Free  air  drag  interference  is  given  in  Figures  32  through  34. 
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Figure  30:  Vectored  Thrust  Lift  Interference 


Figure  33:  Vectored  Thrust,  Drag  Interference,  Vector  Angle  6CP 
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2.3.3  Pitching  Moment  Interference 

Free  air  pitching  moment  interference  also  correlated  well 
with  the  free  air  lift  interference  at  a  given  nacelle  location  and 
nozzle  vector  angle.  This  indicates  that  the  center  of  pressure  of  the 
induced  lift  remains  constant  with  angle  of  attack  for  a  given  nacelle 
configuration.  For  pitching  moments,  the  important  length  parameter  is 
the  distance  from  the  center  of  presure  of  the  induced  lift  to  the  moment 
center.  Therefore,  for  the  pitching  moment  interference  in  free  air,  ' 
Figures  35  through  37  ,  the  nozzle  location  has  been  given  as  the  dis¬ 
tance  from  the  center  of  the  nozzle  exit  to  the  moment  center.  For  a 
swept  wing,  the  average  nozzle  location  is  used. 

2.3.4  Downwash  Interference 

The  effect  of  vectored  thrust  on  downwash  is  shown  in 

Figure  38. 

SAMPLE  PROBLEM  -  VECTORED  THRUST,  FREE  AIR 
c i  *  5.46®  (estimated  power-off  aerodynamic  characteristics) 
a  -  30° 

CL  =  2.4 
C0  -.4062 

Cram  "  0  (model  with  blowing  nozzles) 

Cj  -2.0 

C  -  -  .5356 
m 

c  -  11.179  in. 


Xg  -  -.066  in. 

Z  =  +  2.787  in. 

SL 

x/c  =  . 35 


C. 


=  (-.15+0)  (kr) 


-.15 


With  equation  2.3-1 


-  2.4  -.15  +  2.0  sin  (30  +  5.46) 
*  3.43 

from  wind  tunnel  te3t  data 

-  3.63  @  u  »  5.46 

i*rag 


from  chart  Figure  32  at  CT 


read 


JINT 


-.010 


INT 


calculate  with  equation  2.3.2 

Cjj  *  .4062  -  .010  -  (2.0)  cos  (30  +  5.46)  +  0 
CD  -  -1.2852 

observed  from  TAI  test  data  at  Cj  =  3.43 
CD  .  -1.28 


Pitching  Moment 
from  Figure  35  at  C 


read 


JINT 


C 

m 


+.0450 


INT 


calculate  power  on  with  equation  2.3-3 

C  -  -.5356  +  .0450  +  2.0  sin  30°  +  cos  30°] 

m  11.1/9  11.1/9 

=  -.074 

C  observed  at  C.  =  3.43  wind  tunnel  test 
m  L 


C  =  -.190 
m 

Downwash 

from  Figure  38  read 
Ae  =  +.09° 
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2.4  Vectored  Thrust  in  Ground  Effect 

Vectored  thrust  Interference  effects  in  the  presence  of  the 
ground  were  also  obtained  froa  STAI  wind  tunnel  test  BVWT  099.  Figure 
39  presents  a  comparison  of  free  air  test  data,  free  air  test  data 
corrected  for  ground  influence  and  test  data  in  ground  effect.  These 
show  a  good  correlation  between  the  corrected  data  and  the  test  data 
in  ground  effect. 

2.4.1  Lift  Interference 

As  in  the  case  of  the  power-off  ground  effect  procedure,  the 
C.  vs  ot  curve  in  ground  effect  is  determined  from  the  free  air  curve  by 
adjusting  both  and  ot.  Lift  interference  due  to  vectored  thrust  in 
ground  effect  is  the  sum  of  the  lift  interference  due  to  vectored  thrust 
in  free  air  and  an  additional  increment  for  the  effect  of  ground  proxim¬ 
ity.  This  additional  increment  is  presented  in  Figure  40  .  The  angle 
of  attack  adjustment  is  the  same  as  for  the  power  off  case  (Eq.  2.2-5), 
but  must  be  based  on  C. 

LNET 

Lift  in  ground  effect  with  vectored  thrust  is 


*•  Clge  *■  Cv-visiT  v 

Pouceroff-  f  vecaur 


nCi 

L  IKJT 
around 
e$Te.c^- 
increwvtnl 


+  CT  (<*.  +  «') 

(2.4-1) 


2.4.2  Drag  Interference 

Drag  interference  due  to  vectored  thrust  in  ground  effect  is 
the  sum  of  the  drag  interference  due  to  vectored  thrust  in  free  air  and 
an  additional  increment  for  the  effect  of  ground  proximity.  The  addi¬ 
tional  increment  is  presented  in  Figure  40. 

Drag  in  ground  effect  with  vectored  thrust  is 


c 


D 


=  VN4T 

power  oft  avr 


qroond 

eATecV 

irvcre'nrvenV 


(2.4-2) 
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Thruft  Vactor  Angl#,  <r 


Figure  40:  Change  in  Thrust  Interference  Effects 
Due  to  Ground  Effect 


2.4.3  Pitching  Moment  Interference 

Pitching  moment  interference  due  to  fectored  thrust  in  ground 
effect  is  the  sum  of  the  pitching  moment  interference  due  to  vectored 
thrust  in  free  air  and  an  additional  Increment  for  the  effect  of  ground 
proximity.  This  additional  increment  is  presented  in  Figure  40. 

Pitching  moment  in  ground  effect  with  vectored  thrust  is 

*C»0,k>T  CO*  fl) 

q round 
etCecV- 
mcr«me.riV 

coat*.')  (2.4-3) 


C'rn  - 


C»r\ 

power 


v-np»lK-T 

freeavr 


tn\  c 


"S»vr«jk 


Z*. 
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2.4.4  Downwash  Interference 

Analysis  of  the  test  data  did  not  show  significant  changes  in 
downwash  angle  in  ground  effect  with  the  addition  of  vectored  thrust. 

SAMPLE  PROBLEM,  THRUST  INTERFERENCE  IN  GROUND  EFFECT. 

h/b  -  .208 

a  *  30° 


From  sample  problem  in  Part  2.2  the  test  conditions  in  ground  effect, 
power  off 


V  "  1>93 

GE 

°GE  -  .93 


.282 


-.5088 


The  free  air  vectored  thrust  corrections  at  a  55  30°,  CT  =  0,  nacelle 
x/c  =  .35. 
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C  -  +.065 
miNT 

For  this  example  the  thrust  Interference  effects  In  ground  effect  are 
zero.  Coefficients  in  ground  effect  are  then,  Lift  equation  2.3-1 

CL  -  1.93  -  .2  +  2.0  (sin  30.93) 

-  2.76 

Drag  equation  2.3-2 

Cp  -  .282  -  .025  -  2.0  (cos  30.93)  +  0 

-  -1.458 


Pitching  Moment  equation  2.3-3 

C  -  -.5088  +  .065  +  .255 
m 

=  -.1888 

The  comparable  test  values  at  this  angle  of  attack 

CL  =  2.88 


C  *  -.143 
m 

2.5  Trim 

Any  complete  set  of  longitudinal  data,  lift,  drag,  pitching 
moment,  and  downwash  at  the  tail  may  be  reduced  to  trimmed  lift  and  drag 
by  the  methods  presented  in  this  section.  Mote  that  these  methods  are 
valid  for  relating  long  tail  arms;  close  coupled  tails  or  canards  would 
require  a  considerably  more  involved  analysis. 

2.5,1  Trimmed  Lift 


The  lift  increment  required  to  trim  is  the  increment  required 
at  the  horizontal  tail  1/4  mac  to  reduce  the  pitching  about  the  center 
of  gravity  to  zero. 


trim" 


Cm 

Kfc 


(2.5-1) 
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2.5.2 


Trimmed  Drag 


The  drag  increment  for  trim  is  considered  to  be  made  up  of 
two  components.  First,  the  inclination  of  the  left  vector  since  it  is 
in  the  downwash  of  the  wing.  Second,  the  tail  drag  both  friction  and  the 
tail  drag  due  to  lift. 

A CDtnm=( 

(2.5-2) 


SECTION  III 
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STABILITY  AND  CONTROL  DERIVATIVE  PREDICTION  METHODS 


This  section  includes  methods  for  predicting  vectored  thrust 
effects  on  stability  derivatives,  and  a  sensitivity  study  to  determine 
the  importance  of  each  derivative.  Methods  are  based  on  wind  tunnel 
data  from  Reference  5.  Accuracy  adequate  fc-  preliminary  design  pur¬ 
poses  is  provided.  This  results  in  a  simple,  quick  method. 

Error  charts  and  tables  are  included.  These  should  be  used 
in  conjunction  with  the  sensitivity  study.  The  reader  should  guard 
against  falling  into  the  trap  of  thinking  of  errors  only  in  terms  of 
"percent  error."  Often  it  is  the  increment  of  error  that  is  important. 
For  instance,  in  predicting  the  tail-off  CRg,  an  error  of  200%  would  be 
insignificant  if  the  actual  value  were  only  -.QOQln/deg  1 .  On  the  other 
hand,  if  the  tall-on  CnR  is  .008'Wieg"1 ,  a  15%  error  might  be  quite 
noticeable. 


3.1  Stability  Derivative  Sensitivity  Study 

It  j.s  important  in  the  study  of  an  airplane’s  stability 
characteristics  to  understand  the  consequences  of  errors  in  estimating 
stability  derivatives.  When  the  sensitivity  of  the  dynamic  response 
to  each  parameter  is  known,  effort  to  improve  accuracy  can  be  expended 
on  the  more  important  derivatives. 

Such  a  sensitivity  study  was  performed  for  the  airplane 
shown  in  Figure  41.*  a  nominal  STOL  approach  condition  of  75  knots 
was  selected,  and  stability  derivatives  were  estimated.  The  derivatives, 
together  with  mass  properties  and  reference  dimensions,  are  given  in 
Table  I.  Derivatives  found  to  be  the  more  important  ones  are  listed  in 
Table  II. 


Angle  of  attack  and  sideslip  derivatives  are  based  on  wind 
tunnel  data  from  Reference  5.  Rotary  derivatives  were  predicted  using 
DATCOM  methods. 

Three  degree  of  freedom  equations  of  motion  for  longitudinal 
and  lateral-directional  stability  were  solved,  using  the  nominal 
derivatives.  Then  each  derivative  was  varied  over  a  range  of  +150%, 
except  in  a  few  cases  where  this  would  have  resulted  in  an  unreasonably 
large  increment. 


*This  airplane  is  the  "Baseline  Configuration"  developed  early  in  the 
STAI  program  and  reported  in  detail  in  Appendix  A  of  Volume  I  of  the 
STAI  Series  (Ref.  12). 
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REFERENCE  GEOMETRY 

S  =  1640  FT2 
b  =  114.5  R 
c  =  15.7  R 
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Figure  41:  General  Arrangement  STOL  Tactical  Transport  —  Model  953-801 


TABLE  I 


Si'abi  iity  Derivatives,  Mass  Properties,  and 
Reference  Dimensions  of  Example  Airplane 


All  angles  are  in  radians. 


cx 

a 

«  1.297 

cx- 

a 

=  0 

cx 

q 

=0  Cx  =  -1.77 

u 

CL 

a 

-  7.84 

CL* 

a 

*  -8.18 

CL 

q 

■  7.54  CL  =  .0341 

u 

C 

m 

a 

1! 

» 

4>* 

VO 

Ov 

C 

m* 

a 

-  -6.06 

c 

m 

q 

-  -32.94  C  -  -.453 

m 

u 

\ 

-  -1.415 

CY 

P 

=  .51 

CY 

r 

=  .02 

C 

nP 

=  .27 

C 

n 

=  -.45 

C 

n 

-  -.35 

p 

P 

r 

% 

=  -.191 

C1 

-  -.59 

C1 

-  1.16 

P 

r 

s  - 

1640  ft2 

W  «  133,000  lbs 

• 

IZ2  *  2,62  x  lo6-slug-f t2 

c  = 

15.7  ft 

I  =  1 

XX 

.26  x  106- 

slug-f t2 

IxZ  ■  1.4  x  105-slug-ft2 

t  =*  114.5  ft.  IYy  =  1.46  x  106-slug-ft2 


Angle  of  attack,  a  *  .182 
Thrust  deflection,  a  =  1.13 
Thrust  coefficient,  C  =  1.72 

J 
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Lateral-Directional  Longitudinal 


TABLE  II 


Stability  Derivatives  With  Important  Influence 
On  Airplane  Stability 


Stability 

Derivative 


Major  Influence 

neutral  point  (Note:  when  was  varied,  C  ^  was 

held  constant  so  the  a.c.  was  moving.) 

neutral  point,  short  period  frequency  and  damping 
ratio,  long  period  frequency  and  damping  ratio 

short  period  damping  ratio 
short  period  damping  ratio 

neutral  point,  long  period  frequency  and  damping 


Dutch  roll  frequency,  spiral  stability 

Dutch  roll  damping  ratio,  spiral  stability 

Dutch  roll  frequency,  spiral  stability 

Dutch  roll  damping  ratio,  spiral  stability 

Dutch  roll  damping  ratio,  spiral  stability 

Dutch  roll  frequency  and  damping  ratio, 
spiral  stability 
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Roots  of  the  longitudinal  characteristic  equation  were  plotted 
on  the  s-plane.  Dutch  roll  node  roots  are  also  presented  on  the  s-plane. 
Spiral  node  tine  constants  were  plotted  versus  the  derivative  being 
varied.  These  plots  are  shown,  and  significant  trends  discussed  in  the 
next  sections. 

3.1.1  Longitudinal 

The  influence  of  angle  of  attack,  aerodynamic  lag,  pitch 
damping,  and  speed  derivatives  is  shown  in  Figures  42,  43,  44,  45 
and  46  .  These  Charts  show  that  the  derivatives  critical  to  an  accurate 
determination  of  the  longitudinal  characteristics  are:  ,  C_  ,  Cra  , 

and  CjMu-  “  q 

Sensitivity  of  longitudinal  characteristics  to  variations  of 
the  pitching  moment  due  to  angle  of  attack,  C_  ,  are  shown  in  Figure  42. 
Even  though  is  negative  (the  a.c.  is  more  than  6%  c  aft  of  the  c.g.), 

the  airplane  is  statically  unstable.  (There  is  a  real  root  in  the  right 
half  plane.)  This  is  due  to  the  large  negative  value  of  Cmu.  As 
is  increased  from  its  initial  value,  the  unstable  root  movesto  the  left, 
toward  the  other  real  root,  while  the  complex  root  moves  upwaid.  The 
short  period  frequency  is  increasing  and  the  damping  decreasing.  At 
about  1.5  times  the  initial  C^,  the  previously  unstable  root  goes  to 
the  origin  and  the  airplane  becomes  neutrally  stable.  (The  c.g.  is  at 
the  neutral  point.)  When  C„  is  further  increased,  the  two  real  roots 
couple  and  form  a  long  period  oscillatory  mode,  the  phugoid.  If  Cffla 
were  further  increased,  the  phugoid  mode  may  go  unstable  but  the  air¬ 
plane  would  still  be  statically  stable  (the  neutral  point  would  still 
be  aft  of  the  c.g.) . 

When  Cjj^  is  decreased  the  short  period  frequency  decreases 
and  the  damping  ratio  increases.  The  unstable  root  goes  more  unstable 
and  the  other  real  root  moves  to  the  left.  At  about  .57  times  the 
initial  C^,  the  short  period  mode  becomes  critically  damped.  (The 
short  period  mode  is  now  described  by  real  roots.)  As  Cn,  is  increased 
more,  one  short  period  real  root  moves  to  the  left  while  ?he  other  one 
moves  toward  the  other  stable  real  root.  At  about  .53  times  the  initial 
these  latter  two  roots  couple  and  form  an  oscillatory  mode. 

It  is  necessary  to  know  Cma  accurately  for  reasons  other  than 
longitudinal  dynamics  considerations.  The  aerodynamic  center  should  be 
known  within  about  +1%  MAC  in  order  to  design  the  tail,  locate  the  c.g. 
envelope,  compute  control  surface  deflections  for  trim  and  maneuver, 
etc.  In  this  case,  a  +1%  MAC  error  in  the  aerodynamic  center  location 
corresponds  to  about  a  +15%  error  in  .  Figure  42  shows  that  a  15% 
error  in  C^  will  only  result  in  about  a  5%  error  in  natural  undamped 
frequency  and  a  .05  change  in  damping  ratio. 

Sensitivity  to  lift  curve  slope,  C^,  and  axial  force  due  to 
angle  of  attack,  Cy  ,  is  also  shown  in  Figure  43  .  Varying  Cy  had  no 
noticeable  effect  on  the  unstable  root  and  only  a  smaix  effect  on  the 
others.  A  large  error,  +50%,  In  Cy  ,  should  cause  no  serious  inaccura¬ 
cies.  It  is  hard  to  conceive  of  a  5o%  error  in  Cl  so  this  derivative 

Ot 
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gch  Condition 
74.6  Knots 


Effect  of  Angie  of  Attack  Derivative  ,  C  ,  on  Longitudinal  Dynamic  Stability 


Effect  of  Angle  of  Attack  Derivatives ,  Cx  and  C,  pn  Longitudinal  Dynamic  Stability 


Approach  Condition 
V  =  74.6  Knots 


Figure  44  :  Effect  of  Speed  Derivatives  on  Longitudinal  Dynamic  Stability 


roach  Condition 
74.6  Knot* 
133,000  Lbs 


Figure  46  ;  Effect  of  Aerodynamic  Lag  Derivative*  on  Longitudinal  Dynamic  Stability 


Condition 


Figure  46  ;  Effect  of  Pitch  Damping  Derivatives  on  Longitudinal  Dynamic  Stability 


was  only  varied  +50%.  A  30%  error  in  lift  curve  slope  would  only  affect 
the  undamped  natural  tteqpency  by  about  6%  and  has  a  negligible  effect 
on  the  damping  ratio.  Its  greatest  effect  appears  to  ire  on  the  real 
roots.  As  CI<a  is  reduced,  the  unstable  root  aoves  to  the  left,  couples 
with  the  other  real  root,  and  fores  the  oscillatory  phugoid  node.  With 
a  large  unstable  value  for  C^,  a  20%  or  30%  error  in  &  could  make 
the  difference  between  whether  or  not  the  airplane  was  statically 
stable.  Keep  in  mind  that  C«  was  held  constant  while  Ct  varied,  so 
changing  CLq  also  implies  a  change  in  the  aerodynamic  tenter  location. 

The  influence  of  speed  derivatives  is  shown  in  Figure  44- 
Large  errors  in  and  Cj^  will  cause  no  problem.  However,  Cm  should 
be  accurately  known,  because  large  negative  values  of  cause  the  air¬ 
plane  to  be  statically  unstable  even  though  the  c.g.  is  ahead  of  the 
aerodynamic  center.  has  only  a  small  effect  on  the  short  period 

mode. 

Powered  lift  airplanes  are  likely  to  have  large  values  of 
Cg^.  In  the  trim  condition  a  large  aerodynamic  pitching  moment  is 
required  to  balance  the  thrust  moment.  If  a  speed  change  occurs  these 
two  moments  change  at  different  rates  causing  a  moment  unbalance.  There 
is  another  component,  to  Cm  ,  due  to  thrust  interference  but  this  is 
generally  small  for  -a  vectored  thrust  airplane. 

Effects  of  aerodynamic  lg£  or  the  a  derivatives,  on  longi¬ 
tudinal  dynamics  are  shown  in  Figure  45.  Cl.  has  no  noticeable  effect. 
The  real  roots  are  not  influenced  by  but  ?he  damping  ratio  of  the 
short  period  mode  appears  to  be  sensitive  to  this  term.  As  C™-  is 
increased  the  damping  ratio  increases  and  at  two  times  the  initial  value 
the  short  period  mode  is  critically  damped.  It  would  be  desirable  to 
know  within  40%  in  order  to  know  the  damping  ratio  within  about  10%. 

Sensitivity  to  the  pitch  rate  derivatives  is  shown  in 
Figure  46.  Varying  Cl  had  no  noticeable  effect.  A  ^200%  error  would 
be  negligible.  However^  dynamic  characteristics  are  sensitive  to  Cm  . 

As  is  increased  from  the  initial  value,  the  short  period  damping4* 
ratio **13  increased  without  much  effect  on  undamped  natural  frequency. 

If  Cm  Is  reduced,  undamped  natural  frequency  and  damping  ratio  both 
are  riduced.  The  real  roots  are  only  slightly  affected,  but  if  Cm 
were  increased  still  further  than  shown  in  Figure  46  a  long  perio3 
oscillatory  mode  would  develop. 

3.1.2  Lateral-Directional 

The  influence  of  sideslip,  yaw  rate,  and  roll  rate  deriva¬ 
tives  on  lateral-directional  dynamics  is  shewn  in  Figures  47  through 
52  .  Derivatives  that  must  be  predicted  with  relative  accuracy  are: 

cng»  ^lg»  cnp>  clp»  Cnr»  an<1  clr* 

Sensitivity  to  variations  in  sideslip  derivatives  are  shown 
in  Figures  47  and  48  .  Cy„  has  only  a  small  influence  on  the  Dutch 
roll  mode  and  practically  no  effect  on  the  spiral  mode.  Large  errors 
in  Cy e  would  not  seriously  affect  the  Dutch  roll  characteristics. 

However,  Ci0  and  Cno  strongly  influence 

p  p 


Approach  Condition 
V  *  744  Knots 
W  «  133,000  Lbc. 
a  *  644000. 

7  *  >6  Deg. 
Cj«1.72 


Initial  Values: 

Cv  «- 1.415  (RAD-1) 

0 

Cn  *  .27  (RAD"1) 

e 

Cf  =  -.191  (RAD*1) 


2.6  ^ 


initial 


0  initial 
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both  the  Dutch  roll  and  spiral  nodes.  It  is  desirable  to  know  both  of 
these  derivatives  within  an  increment  of  +.03  -rad-1,  about  10  to  IS 
percent.  It  is  interesting  to  note  that  variations  in  Cn -affect  mainly 
the  Dutch  roll  frequency  while  Ci-changes  affect  primarily  the  Dutch 
roll  damping  ratio.  When  CQ.  is  reduced  to  zero,  the  Dutch  roll  mode 
is  still  stable  and  the  spiral  mode  becomes  stable.  If  Cj_  is  reduced 
to  zero  the  Dutch  roll  mode  remains  stable  but  the  spiral  mode  gets  more 
unstable. 

Accuracy  of  calculations  relating  to  cross  wind  landings  and 
engine-out  conditions  is  directly  related  to  the  quality  of  the  side¬ 
slip  derivatives.  This  should  be  taken  into  account  when  deciding  on 
the  required  accuracy  of  the  derivatives. 

Figures  49  and  50  show  the  effect  of  roll  rate  derivatives. 
Cy  has  no  effect  on  any  of  the  roots  of  the  characteristic  equation 
an S  for  this  purpose  can  be  ignored.  C«  and  effect  both  the  spiral 
and  Dutch  roll  modes  and  should  be  known*3 within  In  increment  of  +.1  - 
rad”1,  or  about  20% .  The  cross  derivative,  Cn  ,  affects  mainly  Dutch 
roll  frequency  and  the  roll  damping  derivative?  ,  affects  mainly  the 
Dutch  roll  damping  ratio.  When  went  to  zero,  ?he  Dutch  roll  damp¬ 
ing  did  too,  even  though  Cn^  and  C?g  both  have  stable  values. 

The  effects  of  yaw  rate  derivatives  are  shown  in  Figures  51 
and  52.  Again  the  side  force  derivative  has  no  effect.  The  cross 
derivative,  Cir,  affects  both  Dutch  roll  damping  ratio  and  frequency. 

The  yaw  damping  derivative  affects  mainly  Dutch  roll  damping.  Cn  and 
both  affect  the  spiral  mode  with  Cir  having  the  greater  influence. 
ReSucing  would  stabilize  the  spiralrmode  while  reducing  the  Dutch 
roll  damping  ratio.  Cnr  and  Cjr  should  be  determined  within  an  incre¬ 
ment  of  +.1  rad-1. 

3.2  Stability  and  Control 

This  section  presents  a  simple  empirical  method  of  predict¬ 
ing  aerodynamic  interference  effects  due  to  vectored  thrust  on  stability 
and  control  derivatives.  The  method  consists  of  applying  a  thrust 
correction  factor  to  the  tail-off  derivative  and  taking  into  account 
the  power  effect  on  the  downwash,  sldewash,  and  dynamic  pressure  at  the 
tail.  It  i3  assumed  that  the  power-off  characteristics  are  known, 
either  estimated  or  from  wind  tunnel  data.  Correction  factors  are  all 
based  on  wind  tunnel  data.  The  wind  tunnel  data  are  presented  in 
Reference  5. 

All  derivatives  and  coefficients  in  this  section  are  net 
values,  that  is;  direct  thrust  forces  are  not  included. 

It  is  appropriate  to  state  here  some  general  observations 
and  opinions  regarding  the  wind  tunnel  yaw  data. 

o  Spanwise  engine  location  has  a  negligible  effect  on  lift 
curve  slope. 
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Double  Amplitude)'1  (Time  to  Damp  to  Half  Amplitude)' 


-1.0  -.5  0  .5  1.0  1.5  2.0 


cy  .  Ca  ,  C#  (RAD-1) 

Yr  r 

Figure  52  ;  Effect  of  Yaw  Rate  Derivatives  on  Spiral  Mode  Stability 
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o  Engine-out  has  little  effect  on  sideslip  derivatives. 

o  Thrust  has  small  effect  on  sideslip  derivatives  except  at  high 

thrust  deflection  or  when  the  nacelles  are  double  podded  inboard. 

o  Angle  of  attack  effects  on  sideslip  derivatives  are  small, 

although  a  little  greater  at  the  higher  wing  sweep  and  in-ground 
effect. 

o  Chordwise  nozzle  positicn  has  a  negligible  effect  on  sideslip 
derivatives. 

o  Ground  effects,  on  sideslip  derivatives,  are  small  except  at 

high  thrust  settings  with  90  degree  thrust  vector  angle.  There 
is  apparently  a  flow  breakdown  at  this  condition. 

o  Thrust  has  negligible  effect  with  flapr  up. 

Thrust  effects  may  be  magnified  by  having  poor  flow  on  the 
model  at  zero  thrust.  With  the  leading  edge  flap  deflected  70  degrees, 
flow  is  stalled  on  the  bottom  of  the  wing,  so  the  trailing  edge  flaps 
are  "seeing"  stalled  air.  Tuft  studies,  in  the  wind  tunnel,  show  that 
the  trailing  edge  flaps  are  in  turbulent  flow  up  to  about  12  degrees 
angle  of  attack.  Also,  the  lift  curve  slope  is  very  high  at  low  angle 
of  attack,  indicating  something  (probably  the  wing  undersurface)  is 
becoming  unstalled  as  angle  of  attack  increases.  All  of  the  yaw  runs 
were  done  at  angles  of  attack  less  than  or  equal  to  12  degrees.  There¬ 
fore,  the  flaps  never  had  "clean"  air  in  any  yaw  run.  The  engines  are 
located  in  this  stalled  air.  They  are  an  energy  source  that  probably 
tends  to  straighten  the  stalled  flow.  This  might  mean  that  the  power 
effects,  presented  here,  are  merely  increments  tending  to  swing  the 
data  back  to  where  the  power-off  data  would  have  been  if  the  bottom  of 
the  wing  had  not  been  stalled. 

3.2.1  Longitudinal  Stability  and  Control 

This  section  presents  a  method  for  estimating  the  aerodyna¬ 
mic  interference  effect  of  engine  thrust  on  longitudinal  stability  and 
control  derivatives.  This  method  has  an  empirical  basis  and  has  been 
derived  from  the  vectored  thrust  blowing  test  (BVWT  099,  Reference  5). 
Methods  for  predicting  lift  and  pitching  moment  are  also  presented  in 
Section  2.  However,  the  methods  presented  here,  although  less  precise, 
are  more  appropriate  for  preliminary  design  purposes  because  they  3re 
faster. 

3. 2. 1.1  Static  Stability  Derivatives 

The  test  pitching  moment  and  lift  curves  are  quite  nonlinear 
with  respect  to  angle  of  attack.  To  obtain  the  results  reported  here, 
slopes  were  measured  at  8°  angle  of  attack,  which  is  representative 
of  takeoff  and  landing  conditions.  The  method  has  been  compared  to 
test  data  at  a*4°  and  8°  and  agreement  is  quite  good  at  both  angles  of 
attack. 
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Figures  53  through  53  show  the  effect  of  nozzle  location, 
vector  angle,  and  Cj  on  tail-off  lift  curve  slope  and  aerodynamic  center. 
These  are  corrections  which  should  be  applied  to  power  off  and  ac  by 
equations  3.1  and  3.2.  a 

Cl«*  Cl~c,«  ♦(%NC*  <31) 

ac  «  acc,.0  ♦  [(^*1/  0.2) 

where: 

ACj^  -  per  degree 

a.c.  -  aerodynamic  center  shift  in  fraction  of  MAC, 
positive  aft 

Subscripts: 

a  -  means  at  constant  spanwise  nacelle  location 

b  -  means  at  varying  spanwise  nacelle  location 

The  nozzle  chordwise  location  is  the  position  of  the  center  of  the 
nozzle  exit  plane  in  percent  of  the  wing  local  chord,  as  shown  in  the 
sketch  below. 


For  ease  of  application  data  are  shown  for  wing  sweeps  of 
0°,  15°,  and  30°.  Power  effects  were  measured  in  the  wind  tunnel  at 
15°  and  30°  only,  the  0°  sweep  is  an  extrapolation  of  these  data. 

Figures  53  and  54  are  for  the  engines  at  27%  and  43.5% 
semispan  locations.  To  account  for  the  effect  of  different  spanwise 
positions  Figure  55  has  been  developed.  Figure  55  shows  the  effect 
of  mean  spanwise  nacelle  position  (average  between  inboard  and  outboard) 
on  a.c.  The  lift  curve  slope  is  not  affected;  however,  inward  move¬ 
ment  of  the  nacelles  has  a  stabilizing  effect  on  a.c.  shift  due  to 
interference. 

The  increments  obtained  from  these  figures  are  compared  to 
the  wind  tunnel  test  results  at  a  *  4°  and  8°  and  for  several  nacelle 
positions,  both  with  single  and  double  pods,  in  Table  III  and  Figures 
56  and  57. 
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NACELLES  AT  27%  and  43.5%  SEMISPAN 


-20  0  20  40  60  80  100 

Nozzle  Location  (%C) 

Figure  54  Effect  of  Vectored  Thrust  on  Aerodynamic  Center 


Mean  Nacelle  Position  (%b/2) 


Figure  SB  Effect  of  Nacelle  Span  wise  Location  on  Aerodynamic  Center 
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TEST  -  PREDICTION  COMPARISON,  ANGLE  OF  ATTACK  DERIVATIVES 
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at  35%  chord 
(nominal) 


Predicted  Change  in  Lift  Curve  Slope 


°  Pi  P4  -  Engines  at  27  and  43.5%  Semispan 
Nozzles  at  0%  Chord 

D  P2  P5  —  Engines  at  27  and  43.5%  Semispan 
Nozzles  at  35%  Chord 

A  P3  p6  ~  Engines  at  27  and  43.5%  Semispan 
Nozzles  at  70%  Chord. 

k  P5  Pg  —  Engines  at  43.5  snd  60%  Semispan 
Nozzles  at  35%  Chord 

0  P12  -  Double  Podded  Engines  at  27% 

Semispan.  Nozzles  at  35%  Chord. 

0  ^13  ~  Double  Podded  Engines  at  43.5% 

Semispan.  Nozzles  at  35%  Chord 


Figure  56  :  Lift  Curve  Slope  Error 
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-Mj  unu.-t.ixLLjL.x-.  Hgt  vj-aue. 


Predicted  a.c.  ShiffOue  to  Thru*t 
I  Ae*.c«  |n»  ^  %  MAC) 


o  PjP*  -  Engines  at  27  and  43£%  Samispen 
Nozzles  at  OK  Chord 

o  P-fs  -  Engines at 27 and 43£% Semispan 
Nozzles  at  3S%  Chord 

A  P3P6  ~  Engines  at  27  and  43.8%  Semispan 
Nozzles  at  70%  Chord 

*  P5P8  -  Engines  at  415  and  60%  Samispen 
Nozzles  at  36%  Chord 


0  P.2  -  Double  Podded  Engines  at  27% 

Samispen  Nozzles  at  35%  Chord 


Figured  :  Aerodynamic  Center  Error 
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The  interference  drag  ten  -.%a  baa  also  been  derived  from 
BVWT  099  (Reference  S)  test  data.  Interference  drag  is  a  function  of 
lift  interference,  vector  angle,  and  nacelle  chordwise  location  and  is 
presented  in  Figures  32  through  34  .  From  these  figures  an  average 


slope  of 


3CL 


INT 


is  obtained.  This  tern,  when  nultiplied  by 


AC.  frcj  Figure  53  gives  the  AC.  tern:  AC_  ■  .3AC. 
a  INT  a  INT  a  INT  a  INT 


The  vectored  thrust  effect  on  horizontal  tail  input  to  lift 
curve  slope  and  aerodynamic  center  is  caused  by  a  change  in  dynamic 
pressure  and  downwash  at  the  tail.  Power-off  tail  effectiveness  should 
be  corrected  for  thrust  effects  by  Equations  3.3  and  3.4 


ac,-“h  *  (aCl«Jc«> 

Aa.cH-(AacM)Cjio  (^J 


(3.3) 


(3.4) 


; - T7T -  is  given  in  Figure  58.  Downwash  is  based  on  tail- 

h-Ulcj.o 

on,  tail-off,  and  tail  control  power  test  data  from  BVWT  099  (Reference 
5) .  The  downwash  shown  is  the  averaged  value  based  on  wing  sweeps  of 
15®  and  30°  and  on  vector  angles  of  30°,  60®  and  90°.  This  shows  good 
agreement  with  downwash  from  wake  rake  data  obtained  in  BVWT  101 
(Reference  5). 

An  attempt  to  measure  power  effects  on  dynamic  pressure  at 
the  tail  proved  unsatisfactory  because  of  wind  tunnel  instrumentation 
problems.  Figure  59  is  presented  instead,  as  a  representative  example 
of  the  effect  of  vectored  thrust.  This  data  was  extracted  from 
horizontal  tail  effectiveness  tests  at  60°  vector  angle. 

Vectored  thrust  has  an  effect  on  the  horizontal  tail  drag. 
However ,  this  is  only  a  email  increment  and  for  preliminary  design 
purposes  may  be  neglected. 
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30 


40 


leg 


n  Horizontal  Tail  Dynamic  Pressure 


3. 2. 1.2 


Derivatives  with  Respect  to  Forward  Speed 


The  speed  derivatives  ,  Cx  ,  and  Cm  are  a  function  of 
both  direct  thrust  and  thrust  interference.  The  uforce  and  moment 
equations  are: 


Ci  *  ' -ACwwettMce  "  CT  0  (3.5) 

Cx  “  ^  rHTCTOE-«.O^CK  V  ^3)  £05.0*4  rf  (3.6) 

C"=  CmC/«»  +AC"We«fehenC6  +  C,(ETcosr.Xt»r 

(3.7) 


Referenced  to  these  equations  the  speed  derivatives  are: 


(3.8) 


C*U*-2C0  +2Cj 


*iAcDj 


«*  CONST 


(3.9) 


C»nu*“2Cj(ZT  cose-  4*  XT  sinc)  -  2Cy 


dC, 


(3.10) 

«*  CONST 


where  CD 


+ 

0 


acd 

Interference 


X 


T 


distance  from  <  .g.  to  thrust  vector  in  fraction  of  MAC, 

positive  fwd. 


■  distance 
positive 


from  c.g. 
down. 


to  thrust  vector  in  fraction  of  MAC, 
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From  the  above  equations,  the  thrust  interference  terms  are 


AC*. 

UMTC*PCMEMCr 

AC*u, 

AC 


■«.  tel 


s CONST 


UINTEftFBftENCC  T  Jax 


(3.11) 

(3.12) 

(3.13) 


oe*  const 


The  terms 


iMJ 


,  ,  and 

can  be  calculated  from  Equations  3.14  through  3.16. 


iC, 


■  J  J0f*CONST 


atasd 


a  -CONST 


J0t=  CONST 


•35fc--.r  *, 

Vc7 

Vc, 

-.I19(C.  +  AC,  ) 

— -  ^vv,l|NT . — — klNT/C  j *2 

x 


(3.14) 

(3.15) 

(3.16) 


where  [CT  +  ACT  ]  is  obtained  from  Figures  29  and  30 

lint  lint  Cj=2 

Since  this  term  varies  with  Cj  by  the  equation: 

Acl  "  [Clint  +  ACLint]  =  [C,NT.ACLJC,.V7  (3,17) 
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The  term  is  obtained  by  differentiating  with  respect  to  C j . 

is  obtained  by  multiplying  -^°  ,  based  on  Figures  30  through  32 

by  .  U&l  is  obtained  by  multiplying  *•  »  based  on 

Figures  33  through  35  ,  by 

3. 2. 1.3  Pitch  Rate  and  Angle  of  Attack  Rate  Derivatives 

No  testing  was  done  to  evaluate  the  effect  of  vectored  thrust 
on  the  wing  body  contribution  to  the  derivatives  Cj,  ,  ,  Cm^,  and  Cz*. 

However,  this  is  expected  to  be  small,  and  existing'WthSds  to  predict01 
the  power  off  wing-body  damping  should  provide  sufficient  accuracy.  The 
horizontal  tail  contribution  to  pitch  rate  damping  derivatives  Cq  and 
Cz  is  influenced  by  engine  thrust  through  the  change  in  dynamic  ^ 
pressure  at  the  tail.  Power  off  Cj,  and  Cz  should  be  obtained  by  exist¬ 
ing  methods  and  the  tail  contributiSn  shoulS  be  corrected  for  thrust 
effects  by  Equations  3.18  and  3.19. 


W*°  TCj%o 


(3.18) 


C*»H  3(CH„)ci.o(tcJ 


(3.19) 


The  horizontal  tail  contribution  to  angle  of  attack  rate 
damping  derivatives  and  C2^  is  a  function  of  both  the  dynai  'c 
pressure  change  and  the  downwash  change  due  to  vectored  thrust.  These 
derivatives  should  be  predicted  by  existing  methods,  with  the  tail 
contribution  corrected  for  thrust  effects  by: 


C|"«h  [%,J 

<  c»*o 


C**h  hWo  jfir  |fCr50 


(3.20) 


(3.21) 
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3. 2.  1.4 


Control  Derivatives 


The  tail  control  derivatives  C^g,  ^x$£*  an<*  are  a^so  a 

function  of  the  dynamic  pressure  at  the  tail.  Power  off  tail  effective¬ 
ness  should  be  predicted  by  existing  methods  such  as  DATCOM,  and  a  power 
correction  applied  by  Equations  3.22  through  3.24. 


(3.22) 

=(C**Jc,«o  (^J 

(3.23) 

(3.24) 

3.2.2  Lateral-Directional  Stability  Derivatives 

This  section  presents  a  simple  empirical  method  of  predicting 
aerodynamic  interference  effects,  due  to  vectored  thrust,  on  lateral- 
directional  stability  derivatives.  Correction  factors  are  all  based 
on  wind  tunnel  data.  The  ^ind  tunnel  data  art  presented  in  Reference  5. 

No  large  error  would  result  in  the  tail-off  sideslip  deriva¬ 
tives  if  thrust  effects  were  Ignored,  The  data  indicate  that  it  is  only 
in  extreme  conditions,  like  90®  thrust  defection  in  ground  effect,  that 
the  thrust  effects  are  large  on  the  more  important  derivatives  Cn„  and 
Clg  :  This  would  probably  only  be  a  transient  condition  and  for  pre¬ 
liminary  design  purposes  might.be  Ignored. 


3. 2. 2.1  Sideslip  Derivatives 


Thrust  effect  on  sideslip  derivatives  can  be  accounted  for 
by  using  the  following  five  correction  factors: 


rV 


(ipV»  1 

to*'C^ic,Jto  * 


Scio¬ 


to 


,%  X 


Cj*o 


where 


8  -  sideslip  angle 

Cy  ~  side  force  coefficient 

C  =  yawing  moment  coefficient 

n 
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Cj  *  rolling  monent  coefficient 

q  *  dynamic  pressure 

o  *  sidewash  angle 

subscripts 

TO  vertical  tall,  denotes  tail-off 
Cj»0  denotes  thrust  is  zero 

Values  for  these  terns  are  presented  in  Figures  60  and  61  .  Side¬ 

slip  derivatives  are  then  conputed  using  Equations  3.25  through  3.27. 


J  ' 


Cj*0 


(3.25) 

A  |n  Mh 

ly*c ,»o/ W 


S s  S*--  +  «v 


"C..O 


’•  ■fe.K- 


7|* 


26) 


.  (3.27) 


where 


S 

b 

a 

SV 

■v 


=  wing  area 
=  wing  span 

*  vertical  tail  lift  curve  slope 
=  vertical  tail  area 

=  distance  from  c.g.  aft  to  vertical  tail  a.c. 

-  distance  from  c.g.  down  to  vertical  tail  a.c. 

“  ratio  of  dynamic  pressure  at  the  tail  to  free  stream 
dynamic  pressure  at  C.  =  0 

*J 
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Thrust  Coefficient  (C j) 


Vectored  Thrust  Effect  Factors  for  Sideslip  Derivatives  Tail  Off 


Nacaties  at  27% 
and  433%  Sami-Span 


The  biggest  correction  factor  is  for  the  sideforce  derivative 
which  is  the  least  important  of  the  three.  See  Section  3.1,  Stability 
Derivative  Sensitivity  Study.  The  more  important  yawing  moment  deriva¬ 
tive  iiaa  no  correction  due  to  thrust.  The  other  important  derivative, 
rolling  moment  due  to  sideslip,  has  a  correction  factor  of  only  1.17 
up  to  a  thrust  deflection  of  60  degrees.  It  can  be  seen  from  the 
derivative  sensitivity  study.  Section  3.1,  that  these  corrections  are 
not  large. 


Sidewash  data  are  shown  in  Figure  62  .  For  this  particular 

model,  thrust  had  no  influence  so  o/i/der  “  1.  However,  it 

•  C**0 

may  be  too  much  of  a  generalization  to  extrapolate  this  result  to  other 
configurations  so  the  term  is  left  in  the  equations.  In  the  absence  of 
additional  data,  assume  no  vectored  thrust  effect  on  sidewash. 


An  attempt  to  measure  power  affects  on  dynamic  pressure  at 
the  vertical  tail  failed  due  to  wind  tunnel  instrumentation  problems. 

It  is  suggested  that  the  values  given  in  Figure  59  ,  for  the  horizontal 
tail,  be  used  until  more  applicable  data  are  available. 

Table  IV  and  Figure  63  show  typical  errors  resulting  from 
the  application  of  the  correction  factors,  presented  in  Figures  60  and 
61  ,  to  the  power-off,  tail-off  data.  While  the  percent  error  is 

sometimes  large,  the  increment  is  usually  small.  These  errors,  when 
viewed  in  conjunction  with  the  derivative  sensitivity  study  presented 
in  Section  3,  are  seen  to  be  small. 

3. 2. 2. 2  Roll  Rate  and  Yaw  Rate  Derivatives 


No  dynamic  testing  was  done  in  the  wind  tunnel  upon  which  to 
base  any  corrections.  Although  the  sideslip  data  obtained  during  the 
wind  tunnel  test  program  is  not  directly  applicable  to  the  yaw  or  roll 
rate  case,  it  does  provide  a  little  insight  upon  which  to  base  an 
opinion  that  the  effect  is  small. 

The  quality  of  the  roll  damping  derivative,  ,  can  be 
improved  by  multiplying  it  by  the  lift  curve  slope  correc?ion  factor, 
as  given  in  Equation  3.28. 


■H$ 


^Cjto 


(3.28) 


This  correction  is  applicable  because  the  roll  damping  is  proportional 
to  the  local  lift  curve  slope  which  should  be  proportional  to  the  3- 
dimensional  lift  curve  slope.  The  tail  contribution  should  be  ignored 
when  computing  C]_  unless  data  on  sidewash  due  to  roll  rate  are  available. 

The  vertical  tail  contribution  to  the  damping  derivatives 
can  be  improved  by  applying  the  dynamic  pressure  ratio  factor.  Equations 
3.29  through  3.33. 
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Symbol 

a(Peg) 

£x 

0<Oeg) 

h/b 
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Engine-Out 

o 
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&&nu& 

30 

oo 

138 

Right  inboard 

A 

8.0 

JS/SISI 0 

30 

M 

138 

Right  Outboard 

□ 

8j0 

2.0 

30 

OO 

140 

None 

0 

8.0 

2.0 

30 

.242 

141 

None 

0 

0 

&/&I&I0 

30 

.242 

142 

Right  Outboard 

A 

8.0 

S/J6/SI0 

30 

242 

142 

Right  Outboard 

A 

8.0 

A/3/0IS 

30 

242 

143 

Right  Inboard 

D 

8 jO 

0 

60 

242 

144 

— 

L 

8D 

2J0 

60 

242 

146 

None 

O. 

6j0 

2.0 

30 

242 

147 

None 

Figure  62 :  Sidewash  at  the  Vertical  Tail 
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/ 

II 

X-Vo 

(3.29) 

<S,* 

vC,»o 

*C,.o 

(3.32) 

c%  = 

CtVvCIM> 

(3.30) 

S' 

c* 

vc,*o 

*c,*o 

(3.33) 

<S  '  (3  31) 


where 

r  is  the  yaw  rate  angle, 

Pt) 

p  is  the  wing  tip  helix  angle,  — 
subscript  V  denotes  vertical  tail  contribution. 

The  power  effect  on  sidewash  due  to  roll  rate  and  yaw  rate 
is  not  accounted  for,  since  there  are  no  data  upon  which  to  base  a 
correction. 


3. 2. 2. 3  Control  Derivatives 

3. 2.2.3. 1  Thrust  Effect  on  Rudder  Power 

The  effect  of  thrust  on  rudder  power  is  shown  in  Figure  64 
Side  force,  due  to  rudder  deflection,  goes  down  with  thrust  at  8  degrees 
angle  of  attack.  At  20  degrees  angle  of  attack,  sideforce  increases 
with  thrust. 

3. 2. 2. 3. 2  Thrust  Effect  on  Aileron  Power 

With  no  aileron  BLC,  thrust  has  little  effect  on  aileron 
power,  see  Figures  65  and  66  .  However,  Figure  67  shows  that  when 
the  ailerons  are  blown,  the  presence  of  thrust  (C  =2.0)  increases  the 
rolling  moment  coefficient  by  about  .01.  Figure  J  65  shows  that  side¬ 
slip  can  have  a  strong  influence  on  the  effect  of  thrust  on  aileron 
power . 

3. 2. 2. 3. 3  Thrust  Effect  on  Spoiler  Power 

Thrust  has  a  strong  influence  on  spoiler  effectiveness. 

This  is  shown  in  Figures  68  and  69  .  Thrust  effect  is  low  at  zero 
angle  of  attack  and  increases  with  angle  of  attack.  This  is  probably 
because  thrust  induces  more  lift  for  the  spoilers  to  operate  on.  See 
CL  effects  on  Figure  53. 
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3.3 


Engine  Out 


This  section  presents  methods  for  calculating  the  pitching 
moment*  rolling  moment*  end  yawing  moment  due  to  engine  failure  for  a 
vectored  thrust  airplane.  These  methods  are  based  on  test  data  from 
EVWT  099  (Reference  5). 

In  order  to  obtain  the  pitching  moment  on  the  airplane  for 
engine  out  conditions*  the  pitching  moment  for  the  all  engine  case  is 
calculated  using  methods  previously  outlined  in  Section  3. 2. 1.1.  From 
this  the  direct  thrust  pitching  moment  is  subtracted  by  the  Equation: 

(Xt *Nr  ♦  ZT  cosrj  (3;34) 

where 

*  Distance  from  moment  center  to  thrust  vector  in  fraction 
of  MAC*  positive  forward 

Zj,  •  Distance  from  moment  center  to  thrust  vector  in  fraction 
of  MAC,  positive  down 

Figures  70  and  71  show  the  effect  of  engine  failure  on 
rolling  moment  and  yawing  moment.  These  data  are  presented  in  the  form 

and  ACq  are  the  lift  and  drag 

changes  due  to  engine  failure  and  include  both  direct  thrust  and  inter¬ 
ference  effects.  The  rolling  and  yawing  moments  are  calculated  by 
Equations  3.35  and  3.36. 

ruled  “  [^^tailed  #,NCot4’r)  +  1 

ENGINE  U  l  ENGINE  nmUNC-'  ^ 


^"FAILED  *  [AClFAUXDC0»^+<r>  +  ^°IMTfAIUCo]  ^3*36^ 
ENGINE  L  ENGINE  ENGINE3 


t  and 

L 


where 


136 


Nacelle  Location  -  %  Semispan 


Nacelle  Location  -  %  Semjspan 
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APPENDIX 


USERS’  MANUAL  FOR  COMPUTER  PROGRAM 


1.  INTRODUCTION 

1.1  Program  Description 

The  program  calculates  the  aerodynamic  characteristics  of  an 
airplane  with  various  types  and  combinations  of  high  lift  devices.  These 
characteristics  include  the  additional  lift,  drag,  and  pitching  rsoments  of 
the  high  lift  devices  and  the  resultant  triamied  totals.  The  high  lift 
devices  include  leading  edge  .flaps  with  and  without  boundary  layer  control 
(BLC),  single  slotted,  double  slotted,  double  slotted /double  hinged, 
and  triple  slotted/double  hinged  trailing  edge  flaps. 

The  program  is  in  FORTRAN  IV  source  code  compatible  with  the 
Control  Data  Corp.  (CDC)  FTN  compiler  on  the  SCOPE  3.3  operating  system 
for  the  CDC  6000/7000  series  computers.  Eight  general  purpose  mathematical 
routines  from  the  BCS  mathematical  library  used  for  interpolation  (INTAB, 
0CTTAB,  TBL,  SEARCH,  TBLU3,  TERP1,  TERP2,  and  TERP3)  are  included  in  FTN 
binary  form  only.  The  deck  also  includes  the  basic  data  tables  developed 
in  the  main  body  of  this  document. 

The  program  requires  a  number  of  tabular  functions  to  be  main¬ 
tained  on  a  permanent  storage  device  file.  The  contents  of  this  file, 
called  the  "permanent  tables",  actually  form  a  substantial  part  of  the 
data  comprising  the  methodology  defined  in  this  report.  The  permanent 
tables  may  be  permanently  updated  or  temporarixy  modified  if  the  user 
wishes  to  change  the  methodology.  Input  procedures  for  such  changes  are 
beyond  the  scope  of  this  report,  but  are  discussed  in  Reference  13. 

1.2  Input  Sequence 

Any  number  of  cases  may  be  done  on  a  single  computer  run. 

Input  data  for  the  whole  batch  begins  with  a  "starter"  card.  Next 
comes  the  first  case's  title  card.  Data  cards  for  that  case  follow. 

Next  comes  the  title  card  for  the  second  case,  and  so  forth. 

1.3  Output  Summary 

A  printout  of  all  permanent  tables  on  the  storage  device  is 
user  selectable.  When  a  table  is  modified,  it  is  printed  for  verifica¬ 
tion,  whether  it  is  a  temporary  or  permanent  modification.  The  title 
for  each  case  is  printed,  followed  by  all  the  input  data  for  that  case. 

The  coefficients  of  lift,  drag,  and  pitching  moment  are 
printed  following  the  input  data  for  each  case  for  free  air,  with 
ground  effect,  and  vectored  thrust.  A  single  run  may  include  several 
thrust  coefficients,  flap  angles  and  wing/ground  distances.  The  out¬ 
puts  for  ground  effect  are  repeated  for  each  wing  height.  This  output 
is  repeated  for  each  flap  angle,  and  the  entire  output  is  repeated  for 
each  thrust  coefficient. 
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2.  INPUT  CARD  PUNCHING  INSTRUCTIONS 

The  STARTER  card  aust  be  punched  with  one-digit  integers  In 
specified  columns.  Title  cards  are  punched  as  72  columns  of  freely 
arranged  alphameric  data.  With  two  exceptions  noted  In  the  discussion 
below,  all  other  data  are  to  be  punched  ir  10-field  seven  digit  format. 
This  format  requires  all  numbers  to  have  a  decimal  point  but  permits 
arbitrary  placement  in  the  seven  digit  field. 

2.1  Input  Card  Order 

"STARTER"  Card 


1st  Case 

TITLE  CARD 

FLAP  TYPE 

Balance  of  Data 

C01 

C02 

C04-C26 

2nd  Case 

TITLE  CARD 

FLAP  TYPE 

Balance  of  Data 

C01 

C02 

C04-C26 

2.2 

Input  Card  Formats 

"STARiER" 

Card 

Column 

Variable 

Description 

7 

IU 

Input/output  unit  for 
permanent  tabular  data* 

**  2  -  Permanent  disk/scratch  disk 
■  4  -  Magnetic  tape/scratch  disk 

14 

NTBL 

Number  of  permanent 
tables  to  be  modified 
(Use  0)* 

21 

IP 

Print  option 

“  0  -  Formal  output  printed 
=  1  -  Tables  and  formal  output  printed 
•  2  -  Tables,  detailed  intermediate  computations, 
and  formal  output  printed 

28 

ICR 

Permanent  tables  on 

card  reader  (Use  0)* 

72-80  May  be  used  for  identification 


*Modif ication  or  updation  of  tables  requires  nonzero  inputs.  This  situa¬ 
tion  is  discussed  in  Reference  13. 


Card  C01 


Column 

Name 

Description 

1-72 

Title 

Card  C02 

Field 

Name 

Description 

1 

AFLAP 

Flap  type  identification. 

1.  Type  1  single  slotted 

2.  Type  2  double  slotted 

3.  Type  3  double  slotted 

4.  Type  4  triple  slotted 

(Refer  to  Figure  6  for  description 

of  the  flap  types.) 

2 

ATTBL 

Number  of  temporary  table  updates. 
0.  To  use  existing  tables  in 
storage. 

Card  C03 

-  Temporary  tables 

(Not  normally  used.) 

Card  C04 

Field 

Name 

Description 

1 

WGROSS 

Gross  wing  area,  3q  ft 

2 

WREF 

Reference  wing  area,  sq  ft 

3 

SPAN 

Wing  span,  ft 

4 

WPERMT 

Wing  semi-perimeter,  ft 

(For  example,  refer  to  sample  problem. 
Page  5.) 

Card  C05 


Field 

Name 

Description 

1 

CPRT.ME 

Extended  wing  chord  length  normal 
to  extended  wing  half  chord  line,  ft 
(Extended  chord  definition  on  Figure  6.) 

2 

CFLAP 

Flap  chord  Types  1  and  2  flap  normal 
to  wing  half  chord  line,  ft 

3 

AQCORD 

Sweep  angle  of  wing  quarter  chord  line, 
deg 

4 

AHCORD 

Sweep  angle  of  wing  half  chord  line, 
deg 
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Card  COS  (Continued) 


Field 

Name 

Description 

5 

AHLAFT 

Sweep  angle  of  aft  flap  hinge  line, 
deg.  This  is  the  hinge  line  for  Types 
1  and  2  flap  or  the  aft  hinge  line  for 
Types  3  and  4. 

Card  C06 

Field 

Name 

Description 

1 

AD  FLAP 

Number  of  flap  deflections,  maximum 
of  4. 

Card  C07 

(For  Flap  Type  1  or  2) 

Field 

Name 

Description 

1 

DFLP, 

First  flap  deflection,  Type  1  or  2 

X 

flap,  deg 

2 

dflp2 

Second  flap  deflection,  deg 

3 

DFLP3 

Third  flap  deflection,  deg 

4 

dflp4 

Fourth  flap  deflection,  deg 

Card  C07 

(For  Flap  Type  3  or  4) 

Field 

Name 

Description 

1 

DFFWj^) 

First  flap  deflection,  Type  3  or  4 

2 

DFAFX) 

flap,  deg 

3 

dffw2) 

Second  flap  deflection. 

4 

dfaf2) 

5 

dffw3) 

Third  flap  deflection. 

6 

dfaf3) 

7 

DFFW4) 

Fourth  flap  deflection. 

8 

dfaf4) 

Card  C08 

(For  Flap  Type  3  or  4) 

Field 

Name 

Description 

1 


Not  used. 
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Card  C08 
Field 
2 

3 

4 

Card  C09 
Field 
1 

2 

Card  CIO 
Field 
1 

2 

3 

4 

5 


(Continued) 
Wane 
CP FLAP 


CAFT 

AHLFWD 


Description 

Forward  flap  chord  (Includes  aft  flap 
rotated  into  forward  flap  chord  plane) 
ft  -  see  igure  6  for  definition. 

Aft  flap  chord  measured  normal  to  wing 
half  chord  line,  ft 

Sweep  angle  of  forward  flap  hinge  line, 
deg 


Name 


Description 


ETEIN  Distance  from  airplane  centerline  to 

inboard  edge  of  trailing  edge  flap, 
semispans 

ETEOUT  Distance  from  airplane  centerline  to 

outboard  edge  of  trailing  edge  flap, 
semispans 


Name 

CLEDGE 

CHORD 


Description 

Leading  edge  flap  chord  measured 
normal  to  wing  quarter  chord  line,  ft 

Wing  chord  normal  to  wing  quarter 
chord  line,  ft 


DLEDGE  Leading  edge  flap  deflection  normal 

to  flap  hinge  line,  deg 

ELEIN  Distance  from  airplane  centerline  to 

inboard  edge  of  leading  edge  flap, 
semispans 

ELEOUT  Distance  from  airplane  centerline  to 

outboard  edge  of  leading  edge  flap, 
semispans 


Card  Cll 


Field 

Name 

Description 

1 

CLAFU 

Lift  curve  slope,  flaps  up,  per 
degree 

2 

AOFLU 

Angle  of  zero  lift,  flaps  up,  deg 

3 

CLHAXU 

Maximum  lift  coefficient,  flaps  up 

4 

ALPHAI 

Initial  angle  of  attack  for  which 
data  is  desired,  deg 

5 

DALPHA 

Increment  in  angle  of  attack  for  which 
data  is  desired,  deg 

6 

CHDLE 

Extended  wing  chord  length  normal  to 
wing  leading  edge  with  trailing  edge 
flap  extended,  ft 

Card  C12 

Field 

Name 

Description 

1 

CLENLE 

Chord  length  of  leading  edge  device 
normal  to  wing  leading  edge,  ft 

2 

CRDNLE 

Wing  chord  length  normal  to  wing 
leading  edge,  ft 

3 

DWTE 

Increase  in  wing  area  due  to  trailing 
edge  flap  extension,  sq  ft 

4 

DWLE 

Increase  in  wing  area  due  to  leading 
edge  flap  extension  (including  only 
the  area  forward  of  trailing  edge  flaps) 
sq  ft 

5 

WPGROS 

Basic  wing  area  inboard  of  outboard 
edge  of  trailing  edge  flap,  sq  ft 

6 

CULE 

Leading  edge  boundary  layer  control 
blowing  momentum  coefficient 

7 

AEDGE 

Leading  edge  flap  type. 

1.  Shaped  leading  edge  devices. 

2.  Conventional  slats. 
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Card  C13 


Field  Name 

1  WGFLA? 

2  SPANLE 

3  CRDDPM 

4  CDPMFU 
Card  C14 

Field  Name 

1 

2  APQCHD 

3  XAC 


4  S2 

5  DSLE 

6  S3 

Card  Cl 5 


Description 

Wing  area  including  leading  and  trailing 
edge  flap  extension  between  the  inboard 
and  outboard  edge  of  the  trailing  edge 
flap,  sq  ft 

Planform  area  of  leading  edge  device, 
sq  ft 

Wing  chord  including  leading  and  trailing 
edge  extension  normal  to  wing  half  chord 
line,  ft 

Minimum  drag  flapsup 


Description 


Not  used,  leave  blank 

Sweep  angle  of  quarter  chord  with  leading 
edge  extended,  deg 

Longitudinal  location  of  aerodynamic 
center  of  basic  trapezoidal  wing,  ft 
(Note,  all  longitudinal  distances  must 
be  from  a  common  reference  point.) 

Wing  area  between  the  inboard  and  out¬ 
board  edges  of  the  leading  edge  flaps, 
sq  ft 

Increase  in  wing  area  from  extension  of 
leading  edge  flaps,  sq  ft 

Wing  area  between  the  inboard  and  out¬ 
board  edges  of  the  trailing  edge  flaps, 
sq  ft 


Field  Name 

1  ALE 


Description 


Location  of  inboard  edge  of  leading  edge 
flaps: 

1.  Start  at  side  of  body. 

2.  Start  outboard  of  side  of  body. 
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Figure  72:  Nomenclature  for  Leading  Edge  Flap 
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Card  C15  (Continued) 


Field 

Name 

Description 

2 

ATE 

Location  of  inboard  edge  of  trailing 
edge  flaps. 

1.  Start  at  side  of  body. 

2.  Start  outboard  of  side  of  body. 

Card  C16  (For  additional  description  of  geometry  on  Cards  C16  and  C17 
refer  to  Figure  72.) 

Field 

Wage 

Description 

1 

CHDYB 

Streamwise  wing  chord  length  at  inboard 
edge  of  leading  edge  flap,  ft 

2 

XB 

Longitudinal  location  of  leading  edge 
of  streamwise  wing  chord  (CHDYB)  at 
Inboard  odge  of  leading  edge  flaos,  ft 
Measured  from  the  coranon  reference 
point. 

3 

CHDPYB 

Streamwise  chord  length  at  Inboard  edge 
of  leading  edge  flaps  with  leading  edge 
extended,  ft 

4 

XPB 

Longitudinal  location  of  leading  edge 
of  streamwise  chord  (CHPDYB)  with  lead¬ 
ing  edge  flaps  extended,  chord  at 
inboard  edge  of  leading  edge  flaps,  ft 

5 

YB 

Spanwise  location  of  inboard  edge  of 
leading  edge  flaps,  ft 

6 

YM 

Spanwise  location  of  leading  edge  flap 
semispan,  ft 

Card  C17 

Field 

Name 

Description 

1 

CHDYM 

Streamwise  chord  at  mid  span  of  leading 
edge  flap,  ft 

2 

XM 

Longitudinal  location  of  leading  edge 
of  chord  CHDYM,  ft 

3 

CHDPYM 

Extended  streamwise  chord  at  midspan  of 
leading  edge  flap,  ft 

4 

XPM 

Longitudinal  location  of  leading  edge 
of  chord  CHDPYM,  ft 

147 


Card  C18  (For  additional  description  of  geometry  on  Cards  C18  and  C19, 
refer  to  Figure  73.) 


Field 

Name 

Description 

1 

CHDYBT 

Streamwise  chord  at  inboard  edge  of 
trailing  edge  flap,  ft 

2 

XBTE 

Longitudinal  location  of  leading  edge  of 
chord  CHDYBT,  ft 

3 

CDPYBT 

Streamwise  chord  length  at  inboard  edge 
of  trailing  edge  flaps  with  trailing 
edge  extended,  ft 

4 

XPBTL 

Longitudinal  location  of  leading  edge 
of  chord  CHPYBT,  same  as  XBTE,  ft 

5 

YBTE 

Spanwlse  location  of  chord  CHDYBT,  ft 

6 

Card  C19 

YMTE 

Spanwise  location  of  r'dspan  of  trailing 
edge  flap,  ft 

Field 

Name 

Description 

1 

CHDYMT 

Streamwise  chord  at  trailing  edge  flap 
semlspan,  ft 

2 

XMTE 

Longitudinal  location  of  leading  edge 
of  CHDYMT,  ft 

3 

CPYMTE 

Streamwise  chord  at  trailing  edge  flap 
semispan,  includes  trailing  edge  flap 
extension,  ft 

4 

XPMTE 

Longitudinal  location  of  leading  edge 

of  CPYMTE,  same  as  XMTE,  ft 

Card  C20 


Field 

Name 

Description 

1 

CREF 

Reference  chord  length  for  pitching 
moment,  ft 

2 

CFLIB 

Streamwise  flap  chord  at  inboard  edge 
of  trailing  edge  flap,  flap  type  1  or 
2,  ft 

3 

CFLOB 

Streamwise  flap  chord  at  outboard  edge 
of  trailing  edge  flap,  flap  type  1  or 
?,  ft 
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Card  C20  (Continued) 

Field 

Name 

Description 

4 

CPMI3 

Streamwlse  wing  chord  at  inboard 
edge  of  trailing  edge  flap  (includ¬ 
ing  trailing  edge  flap  extension),  ft 

5 

CPMOB 

Streamwlse  wing  chord  at  outboard  edge 
of  trailing  edge  flap  (Including 
trailing  edge  flap  extension),  ft 

Card  C21 

Field 

Name 

Description 

1 

CPFLIB 

Streamwlse  flap  chord  at  inboard  edge 
of  trailing  edge  flap,  includes  aft 
flap  rotated  about  hinge  line  into 
forward  flap  chord  plane,  flap  type  3 
or  4,  ft 

2 

CPFLOB 

Streamwlse  flap  chord  at  outboard  edge 
of  trailing  edge  flap  (includes  aft 
flap  rotated  about  hinge  line  into  for¬ 
ward  flap  chord  plane),  flap  type  3  or 
4,  ft 

3 

XIB 

Longitudinal  intersection  of  CPMIB  and 
wing  half  chord  line,  ft 

4 

XOB 

Longitudinal  intersection  of  CPMOB  and 
wing  half  chord  line,  ft 

5 

XPQCRD 

Longitudinal  location  of  wing  quarter 
mac  with  the  leading  edge  flap  extendei 
ft 

6 

XCG 

Longitudinal  location  of  center  of 
gravity,  ft 

Card  C22 

Field 

Name 

Description 

1 

CMOLFU 

Pitching  moment  coefficient  at  zero 
lift,  flaps  up 

2 

EPTEIB 

Distance  from  airplane  centerline  to 
intersection  of  CPMIB  and  wing  half 
chord  line,  9emispans 
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Card  C22  (Continued) 


Field 

Name 

Description 

3 

Card  C23 

EPTEOB 

Distance  from  airplane  centerline 
to  intersection  of  CPMOB  and  wing 
half  chord  line,  semispans 

Field 

Name 

Description 

1 

THQM 

Height  of  horizontal  tail  quarter 
mac  above  wing  chord  plane,  ft 

2 

TLQH 

Tail  length,  from  wing  quarter  mac 
to  tail  quarter  mac,  ft 

3 

NWINGH 

Number  of  ground  heights  to  be  con¬ 
sidered,  maximum  of  4  (must  be  right 
adjusted  integer  in  Column  21) 

4 

WHQMj 

First  ground  height,  ground  to  wing 
quarter  mac,  ft 

5 

whqm2 

Second  ground  height,  ft 

6 

whqm3 

Third  ground  height,  ft 

7 

Card  C24 

whqm4 

Fourth  ground  height,  ft 

Field 

Name 

Description 

1 

STAIL 

Horizontal  tail  area,  sq  ft 

2 

EPSO 

Downwash  angle  at  horizontal  tall, 
flaps  up,  zero  angle  of  attack,  deg 

3 

DEPDAL 

Rate  of  change  of  downwash  with  angli 
of  attack,  flaps  up,  degrees/degree 

4 

DCPDM 

Horizontal  tail  minimum  drag  coef¬ 
ficient,  referred  to  horizontal  tail 
area 

5 

DCDDCL 

(dC])/dCi,2)j  tail  Induced  drag  factor 
referred  to  tail  lift  coefficient 
and  tail  area 
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Card  C25 


Field 

Name 

Description 

1 

CDRAH 

Ram  drag  coefficient 

2 

ENGVEC 

Engine  vector  angle  referred  to  wing 
chord  plane,  deg 

3 

XNACEL 

Longitudinal  distance  from  wing  lead¬ 
ing  edge  to  nozzle  exit  divided  by 
wing  chord  length  at  the  same  stream- 
wise  station 

4 

XNOZLE 

Longitudinal  distance  from  center  of 
gravity  to  nozzle  exit  centerline,  ft 

5 

ZNOZLE 

Vertical  distance  from  center  of 
gravity  to  nozzle  exit  centerline,  ft 

6 

ZINLET 

Longitudinal  distance  from  center  of 
gravity  to  centerline  of  inlet  face,  ft 

7 

ZINLET 

Vertical  distance  from  center  of  gravity 
to  centerline  of  inlet  face,  ft 

Notei  For  multi-engine  configurations  use  average  values  for  nozzle 
and  inlet  locations. 

Card  C26 

Column 

Name 

Description 

7 

NCG 

Number  of  engine  thrust  coefficients; 
(integer  5)l 

Field 

Name 

Description 

2-NCJ+l 

cj1 

Array  of  engine  thrust  coefficients; 

NCJ  fields  are  used. 


3.0  EXECUTING  FROM  CARDS 

To  execute  the  program  from  the  FORTRAN  source  deck  with  the 
tables  on  cards,  saving  the  tables  on  a  magnetic  tape  file,  TAPE4,  and 
printing  the  tables,  the  deck  setup  required  is  shown  on  Figure  74. 

3.1  Using  the  External  Table  File 

The  deck  is  stacked  with  the  REQUEST  control  card  as  in  Section 
3.0  for  magnetic  tape.  The  tables  are  retrieved  from  the  tape  without 
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Figure  74:  Program  Deck  Stacking 


L  53 


the  full  set  of  tables  in  the  data.  The  Installation  convention  for 
identification  of  the  tape  is  required  on  the  REQUEST  card.  A  peruar.ent 
disk  file  can  be  used  Instead  of  the  magnetic  tape  file  and  peraanent  or 
tenporary  table  updates  nay  be  supplied  as  discussed  in  Reference  13, 
Section  2.  Coluan  28  of  Card  1  is  zero. 

If  execution  without  the  FORTRAN  source  deck  is  desired,  con¬ 
sult  a  programmer  for  the  most  desirable  approach  for  your  installation. 
This  would  be  preferred  ina  production  situation  where  no  changes  are 
being  made  to  the  program. 

4.  SAMPLE  INPUT  AND  OUTPUT 


The  sample  problem  used  in  this  document  is  presented  on  an 
input  data  form  in  Table  V,  followed  by  the  corresponding  output. 

The  output  data  begins  with  a  recapitulation  of  the  input  data. 
The  constants  defining  the  power  off  lift,  moment  and  drag  curves  are 
then  stated,  followed  by  tabulated  tall-off  lift,  moment  and  drag,  and 
by  tabulated  trimmed  lift  and  drag,  for  both  power-off  and  power-on 
conditions.  If  ground  effect  data  is  called  for,  tabulated  character¬ 
istics  in  ground  effect  come  next. 


5.  PROGRAM  LISTING 

A  FORTRAN  IV  source  code  listing  of  the  program  begins  on 

page  163. 
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TABLE  V 
SAMPLE  INPUT 


TRIPLE  SLOTTED  FLAP  TYPE  4  SAMPLE  PROBLEM 
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